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Abstract 



The radiative correction to the Lamb shift of order o^iZoifm induced by the 
hght by hght scattering insertion in external photons is obtained. The new 

contribution turns out to be equal to — 0.122(2)Q;^(ZQ;)^/(7rn^)(mr/m)^m. 
Combining this contribution with our previous results we obtain the complete 
correction of order cp-(Za)^m induced by all diagrams with closed electron 
loops. This correction is 37.3(1) kHz and 4.67(1) kHz for the \S- and 25*- 
states in hydrogen, respectively. 



1 INTRODUCTION 



Recently we started a calculation of all contributions to the Lamb shift of 
order a'^{Za)^m. It was shown that there exist six gauge invariant sets of 
diagrams, which produce such corrections. All these diagrams may be ob- 
tained by different dressings from the skeleton diagram which contains two 
external photons attached to the electron line. Contributions induced by po- 
larization operator insertions in external photons, by simultaneous insertion 
of a radiative photon in electron line and a one-loop polarization operator 
in the external photons and by polarization operator insertions in radiative 
photons have been calculated in 0, ||, . 

We present below a calculation of the a'^{Za)^m contribution to the Lamb 
shift in hydrogenlike ions, induced by penultimate gauge invariant set of 
graphs. These are the graphs containing one-loop light by light scattering 
insertions between two external photon lines (see Fig.l). 

The contribution to the Lamb shift, produced by the diagrams in Fig.l 
is given by the expression (see, e.g. lH) 

. ^ a'^(Za)^ 1 [ d^q 1 f d^\i , . „ ^^ 

irn-^ m tt'^ J tc^i J 47rk* 

where is the four-momentum of the upper photon lines, = (0, k) 
is the spatial momentum of external photons, Aaf3 and Ba/s are electron-line 
factors, and Sa/soo is the light by light scattering tensorQ. 

2 SIMPLIFICATION OF THE EXPRESSION 
FOR THE ENERGY SHIFT 

We begin the calculation with a simplification of the integrand in eq.(|ip. The 
general expression for the light by light scattering tensor is 

d'^p \ _ f d'^p 



SaPQO — I -^{"^La/BOO + Capoo) = [ —^SapOO-, (2) 



where 



^Really Sapm is defined as one fourth of the hght by light scattering tensor; see nor- 
mahzation conventions for this tensor as weU as for the electron factor in the next section. 
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corresponds to each of the ladder diagrams a and b in Fig.l, and 

Cq/300 = (4) 

Tr{7o(]3 - k + m)^p{p - k - q + m)^Q{p - q + m)^a{p + fn)} 



4:DiD2D3Di 



corresponds to the crossed diagram c in Fig.l. Factors in the denomina- 
tors above are defined as follows 

Di = p^ — w? , D2 = {p~kY — m^ , D-i = {p — qY — rn?, = {p — k — qY — rn?. 

(5) 

The explicit expression for the electron factor in eq. (|I|) is the sum of two 
terms, where 

. Tr{(7o + l)7a(7o^ + 9 + ^7/3} 2 (a\ 



g2 + 2mgo q"^ + 2mgo 

corresponds to the diagram with nonintersecting upper photon lines and 

Tr{(7o + l)7/3(7om-g + m)7a} 2 /n\ 
" 4(g2 - 2mqo) ^ 



g2 — 2mgo — 2mgo 

corresponds to the diagram with crossed upper photon lines 0. 

Note that Aa/siq) = Bap{—q) and both Aap and Bap are symmetric rel- 
ative to permutation of indices. Hence, only the even in q and symmetric 
under permutation of indices a and (3 part of the light by light scattering 
tensor is relevant for calculations below. Due to gauge invariance terms in 
the light by light scattering tensor which are proportional to vectors qa and 

^We have included an additional factor of m? in the definition of the electron line factor 
simply to preserve the appearance of the expression for the energy shift in eq. (|l|) even after 
transition to the dimensionless momenta. 
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are also irrelevant for calculation of the contribution to the Lamb shift 
and we omit such terms below. 

Since Aa/siq) = Bap^—q) the expression for the energy shift may be pre- 
sented in the following form 



= ^ J m — f- / / A^p[q)[S^p,,[q) + S^poo{-q)] 

irn-^ m it'' J irH q^ J 47rk* 

(8) 



a'^(Za)^ ^rrir^.. 1 r d^q 1 r (Pk 

mi — Y— / — — - / — — X 



7m^ m vr^ J ir'^i q^ J 47rk^ 

-Y^Aaf3{q)[Sai3oo{q,k,p) + So,poo{-q, k, p)], 

where 5a/3oo is the integrand of the light by light scattering tensor sym- 
metrized in a and j3. 

Both k and p are dummy integration momenta and we use substitutions 
k —k and p — > —p while calculating the second product in the integrand in 
eq.(P). After these substitutions denominators in the integral representation 
for the light by light scattering tensor with external momenta of opposite 
sign return to the original form (see eq.(^) and eq.(^)). Moreover, it is easy 
to see that 

SapQQ{-q,-k,-p;m) = SapoQ{q,k,p; -m), (9) 

if one recollects that there exists an additional argument m in the light by 
light scattering tensor and takes into account explicit expressions in eq.(|^) 
and eq.(^). Traces in the numerators on the right hand sides in these equa- 
tions convert these numerators into even polynomials of m and, hence, the 
expression for the energy sphtting in eq.(^) reduces to 

/^E = (10) 

a'^iZaf ,mr,.2 r I rd^k r d^p 
-m' ' ' ' 



TTV? m Tc"^ J Ti'^i q^ J Ank'^ 



/dp 
-^Aa/siq) Saf^oo {q, k,p;m). 



CALCULATION OF THE LIGHT BY LIGHT 
SCATTERING TENSOR 
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3.1 LADDER DIAGRAM 

The explicit expression for the tensor in eq.(|^) symmetrized over indices a 
and /?, has the form 

= - n2n n {-Mp'^9ap + 2{kp){pq)g^p + A{kp)p„pp (11) 

+ {kp)gapm^ - {kq)p^gap + {kq)gapw? - k^p^pp + kaP^m^ - kpp^p^ + k^Pain^ 

+P'^9aP -P^ iPq)ga/3 - '2p^plgaf3 " 2p^PoPag/30 " '^^P^PoPpQaQ - ^p^PaQogafS " ^p^PaPp 

-2p^gapm'^ + 4:{pq)plgaf3 + {pq)gai3'm^ + ^pIpo^Pp + 2plga(^m^ + 2pop 

+2poPi3gaom^ + 2poqogc,pm^ + 2paPi3m^ + gapm^- 

Next we combine denominators with the help of the Feynman parameters 
X and u 

(1 - u)Di + u[xD2 + (1 - x)D3] = _ A, (12) 
p' = p — kux — qu{l — x), 
A = —k'^xu{l — ux) — q^u{l — x)[l — u{\ — x)] + 2{kq)u^x{l — x) + m^, 

Thus we obtain 

^0/300 = (13) 



-6 ^ ^ rf«n(l -u)J Jr{-P"L% + -p'^nX + ^Soo) 



where 



-^iflOO — ~9a0gi30 + (14) 



-^Soo = -2Pgaogf3o{uxf + APg^piuxf - k'^g^piux) (15) 
A{kq)gaogpQ{ux)[u{l - x)] + %{kq)gc,p{ux)[u{l - x)] - {kq)gc,p{ux) 
-{kq)gai3[u{l - x)] - {kq)gcp - Ak^kp{uxf - 2q^gc,QgiiQ[u{l - x)f 
+Vfi'a/3[^(1 - x)f - q^ga(5[u{l - x)] - Aqlgo,i3[u{l - x)f 



-'2Qo9ai3[u{l - x)] + 2gaQgpQm^ - 2gai3m'^, 

-^Soo = k'^QccpiuxY - k'^Qapiuxf + Ak"^ {kq)gai3{uxf[u{l - x)] (16) 

-k'^{kq)gap{uxf - ?,k'^{kq)gaf5{uxY[u{l - x)] + k'^{kq)gap{uxf 
—2k'^kaki3{uxY + 2k'^kakp{ux)^ — 2k^kaqogi3o{ux)^[u{l — x)] 
-2k'^kpqogaoiuxf[u{l-x)]+2k'^q^ga/3{ux)^[u{l-x)]^-k^q^ga/3iux)^[u{l-x)] 
-k'^q^ gaf:i{ux)[u{l - x)]^ + 2k^q^ gap{ux)[u{l - x)] - 2k^qQgap{uxf[u{l - x)]"^ 
—2k'^qQgai3{uxY[u{l — x)] — 2k'^ ga^m? {uxY + k'^gap'^n^iux) 
+A{kqfg^p{uxf[u{l - x)f - 2{kqfg^p{uxf[u{l - x)] 
-2{kqfgai3{ux)[u{l - x)f - A{kq)kc,k/3iuxf[u{l - x)] 
-4:{kq)kaqogpoiuxY[u{l - x)]^ - 4:{kq)kpqogao{uxf[u{l - x)]^ 
^A{kq)q^gap{ux) [u{l-x)f -?,{kq)q^gap{ux)[u{l-x)f - {kq)q^gap[u{l-x)f 
+ (^?)?^5'a/3['ii(l - a;)]^ - 4(/i;g)go5r„^(iix) - x)]^ - 4(/i;g)5r„;jm^('ua;) [^(1 - x)] 

+ {kq)gapm'^{ux) + {kq)gapm'^[u{l - x)] + {kq)gapm'^ 
—2ka.kj3q^{uxY[u{ \ — x)]^ — 2kakpq^{ux) \u{\ — x)]^ + SkakpqQ{uxY[u{l — x)]^ 
+2kaki3W?{ux)'^ + 2kakjsm^{ux) — 2kaq^qQgpQ{ux)[u{l — x)]^ 
+2kaqQgi3Qm'^{ux)[u{l - x)] - 2k/3q'^qogao{ux)[u{l - x)f 
+2kfsqogaom^{ux)[u{l - x)] + q^^g^pluil - x)]^ - q'^gafs[u{l - x)]^ 
-2q^q^o9c.M^ - ^)]' + 2g'goV/3[M(l - x)f - 2q^g^pm'[u{l - x)]^ 
+q^gai3m^[u{l - x)] + 2qQg„3m^[u{l - x)]^ + 2qQgapm^[u{l - x)] + ga/sm'^- 

Next we perform the Wick rotation and integration over the loop momen- 
tum p'. The momentum integral diverges logarithmically. It is well known 
that even this log of the ultraviolet cutoff cancels if one adds the crossed 
diagram contribution. A finite term breaking gauge invariance remains as 
the only remnant of the divergence. Since it is still necessary to subtract the 
value of the light by light scattering tensor at zero momentum in order to 
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restore gauge invariance we will now perform the subtraction directly in the 
ladder diagram contribution. 

Integration is performed trivially and we obtain 

^ a/300 = (17) 

Next we do the subtraction. It is convenient first to perform auxiliary 
integration by parts. This gives one a chance to separate a subtraction term 
connected with the logarithm and to get rid of this log in the integrand 
simultaneously. We use the identity 

n M (1-^)2(2^ + 1) , A2 
duu{l-u) log— = ^ ^ 

- ^(1 - uy{2u + l){-ex{l - 2xu) 
Jo 6A 

-g^(l - a;)[l - 2u{l - x)] + 4:{kq)ux{l - x)} 

= ^ log ^ - /' ^(1 - ufi2u + l){-exil - 2xu) 
JO oA 

-g^(l - a;)[l - 2u{l - x)] + 4:{kq)ux{l - x)} 



to deal with the log on the right hand side in eg. (p!7|) . 

Subtraction of the logarithmic term at zero momenta is performed simply 
by omitting the log as well as the constant term after integration by parts. 

Some finite terms in eq.([T7|) also need subtraction. As was mentioned 
above it is necessary to subtract the value of the ladder contribution to the 
light by light scattering tensor at zero external momentum. Considering 

explicit expressions for the functions L^^poo eq. ([T^) it is easy to see that 

only terms 2gaogi3om? — 2gaf3m'^ in function L^^loo ^"^^ term gaprn^ in 

function -Lq,^oo need subtraction. Now, subtraction is performed simply by 
the replacements 

1 ^ -1 ^ (19) 

A A(A;,g) A(0,0) ^ ' 

k'^xuil — ux) + ^^^(1 — x)[l — u{l — x)] — 2{kq)u^x{l — x) 
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A2^A2(fc,g) A2(0,0) ~ ^A(fc,g) A(0,0)^^A(fc,g) ^ A(0,0)^ 

= (2m^ — k'^xuil — ux) — q^u{l — a;)[l — u{l — x)] + 2(kq)u^x{l — x)) 

k'^xu{l — ux) + q'^u{l — x)[l — u{l — x)] — 2{kq)u'^x{l — x) 

m^A^ 

in the denominators of the terms just discussed. 

We obtain then the expression for the subtracted ladder diagram contri- 
bution 

>CSo = (20) 
Jo Jo 3 A A2 

where 

-^S? = -2k'^gaog0ouvx{l - u) - 6k^gaogpouv{ux) (21) 
+2k'^gaogi3o{uxf{i - u) + 6k'^gaogpou{uxf - k'^gaogpovx{l - u) 
+k'^gaog/3ox{ux){l -u) + 2k^gai3uvx{l -u) + 6k'^gai3uv{ux) 
-2k'^gai3{uxY{l -u) - l2k'^gai3u{uxY + 3k'^gapu{ux) + k'^gapvx{l - u) 
-k^gai3x{ux){l -u) + S{kq)gaogpo[u{l - x)]{ux){l - u) 
+24:{kq)gaogf3ou{ux)[u{l - x)] + 4:{kq)gaogf3o{l - x){ux){l - u) 
—'&{kq)gai3[u{l — x)\{ux){\ — u) — 'iQ{kq)gapu{ux)[u{l — x)] 

+3{kq)ga(}u{ux) + 3{kq)gaf3u[u{l - x)] + 3{kq)gapu 
—4:{kq)gap{'i- — x){ux){l — u) + 12kaki3u{uxY + 6kaki3u{ux) 
-2g^5'ao5'/3os[«(l - -u)- 6q'^gaog(3oSu[u{l - x)] 
-q'^gaog (303(1 - x){l - u) + 2q^gaog/3o[uil - x)][u{l - - u) 

+Qq'^gaog/3ou[u{i -x)]'^ + q^gaog^oi.^ - x)[u{i - x)]{i - u) 

+2q^gaps[u{l - x)](l - u) + e>q^ gai3su[u{l - x)] + q^gaps{l - x){l - u) 
-2g V/3[m(1 - x)] [u{l - x)] {l-u)- I2q^g^pu[u{l - x)f + 3q^gapu[u{l - x)] 
-q'^gafsil - x)[u{l - x)]{l -u) + Uq^g^puluil - x)f + 6qQga(3u[u{l - x)], 
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^aloo = -k'^Qapv^iuxf + k^Qapiuxf - k^gaf^{uxf (22) 
+4:k'^{kq)gapv{uxf[u{l - x)] + 4:k'^{kq)gap{uxf[u{l - x)] - k'^{kq)gap{uxf 
—3k'^{kq)gai3{ux)^[u{l — x)] + k'^{kq)gap{ux)'^ — 2k'^kakp{uxY + '2,k'^kakp{ux)^ 

-2k'^kaqogf5o{uxf[u{l - x)] - 2k'^kpqogaoiuxf[u{l - x)] 
—2k'^ q"^ ga/3sv{ux)[u{l — x)] + 2k'^q'^ ga/3{uxY[u{l — x)]"^ — k'^ q"^ gai3{uxY[u{l — x)] 
-k'^q'^ga(}{ux)[u{l - x)]'^ + 2k'^q'^g^p{ux)[u{l - x)] - 2k'^qlg^f){uxf[u{l - x)]'^ 
-2k'^qQga/3{ux)'^[u{l - x)] + 2k'^gapm'^v{ux) - 2k^gai3m^{uxf + k'^gapm'^{ux) 
-2{kqfgai3{uxy[u{l - x)] - 2{kqygc,p{ux)[u{l - x)]'^ 
-4{kq)kaki3{ux)^[u{l - x)] - 4:{kq)kaqogpo{uxy[u{l - x)f 
-A{kq)kf)qogao{uxf[u{l - x)]^ + A{kq)q'^gafss{ux)[u{l - x)f 
+A{kq)q^gai3iux) [u{l -x)]^- 3{kq)q^gai3iux) [u{l- x)]'^ - {kq)q^ga0[u{l- x)]^ 
+ {kq)q'^gaf3[u{l-x)]'^-A{kq)q'^gapiux)[u{l-x)]^-8{kq)gapm'^iux)[u{l-x)] 
+ {kq)gapm'^{ux) + {kq)gapm^[u{l-x)] + {kq)ga(iW? -2kakpq^{uxf[u{l-x)f 
—2kaki3(f{ux)\u{l — x)^ + %kaki3q^{uxY[u{l — x)^ + 2kaki3m^{uxY 
+2kakf3m^{ux) - 2kaq^ q^g pq{ux)[u{1 - x)f + 2kaqQgpQm^{ux)[u{l - x)] 
-2ki3q^qogao{ux)[u{l - x)f + 2kpqogaom^{ux)[u{l - x)] - q^gaps^[u{l - x)f 
+q'^gap[u{l - x)]^ - q'^gapiuil - x)]^ - 2q'^qQgap[u{l - x)]^ 
+'^q\o9a(3Hi - x)]^ + 2q^g^pm^s[u{l - x)] - 2q^g^pm^[u{l - x)]^ 
+q'^gai3m'^[u{l - x)] + 2qlga0m'^[u{l - x)f + 2qlgaf3m'^[u{l - x)], 
V = 1 — ux, s = 1 — u{l — x). 
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3.2 CROSSED DIAGRAM 

The symmetrized expression for the crossed diagram tensor in eq. @ has the 
form 

Capm = - n n n {^^ Vfl'aO^/jo - k'^P^Oaf^ - 2k'^{pq)gaogf3o (23) 

+ k^{pq)gal3 - '^k'^PoPagpO - '^k'^PoPpgaO + '^k^PaP/S + k^PaQog^O + k^P^QogaO 

-2P gaogpom^+k^ ga^m^ -'2{kp)kaPp-2{kp)ki3Pa-4:{kp)p^ gaogpo+'2{kp)p^ Qap 
+4:{kp){pq)gaogf3o-2{kp){pq)gai3+4:{kp)poPagi3o+4:{kp)poPf^gao-2{kp)p 
-2{kp)pf}qogao - 2{kp)q'^gaog^o + {kp)q^gap - 2{kp)qlgaf} + 4:{kp)g 
-2{kp)gaf3m'^ + ikq)kaPp+ikq)kf3Pa+2ikq)p'^gaogpo-{kq)p'^ga(i-2ikq)poPag(io 

-2{kq)poPpgao + 2{kq)poqogap - 2{kq)gaogpQm^ + {kq)gapm^ + 2kakpp^ 
-2kakf3{pq) - 2kakpm^ + 2kaP^pQgpo - kap'^qogpo - 2ka{pq)pogpo - "^ko.plpp 
+4:kaPoP(iqo + kaPoq'^gpo - 2A;„po5'/3o"^^ + kaqogfiom'^ + 2A;^pVofl'aO - ^/3P^go5'aO 
-2kp{pq)pogao - 4:kpplpa + '^kppQp^qo + k^p^q^g^o - 2kppQgaQm^ + kpqogaom^ 

+2p'^gaOgf30 - P^gaP - 4p^(pg)5'aOfi'/30 + 2p'^{pq)gaf3 - 4p^PoPa5'/30 " ^P^PoPpg^ 

+2p^Pa?o5'/30 + 2p'^Pf3qogaO + 2p'^q'^gaogi30 - P^q^gap + 2p^qlg^p - ^p^g^gpom^ 

+2p^go,pm^ + 4:{pq)poPagi3o + 4:{pq)poPi3gao - '^{pq)Poqogaf3 + '^{pq)g 

-2{pq)gapm'^ + SplpaPp + 2plq^gap - SpoPaPpqo - 2poPag^fi'/3o + 'ipoPagpo'm'^ 

-2poppq'^gao + Ap^ppg^Qm^ - 2pc.gofl'/3o"^^ - 2ppqogao'm'^ - 2g^5'ao5'/3o"^^ 

+q^gapm^ - 2qlgapw? + 2gaQgpom^ - gapn/]. 

Next we combine denominators and perform momentum integration. Un- 
hke the case of the ladder diagram there are four different denominators now 
and an extra Feynman parameter z is needed to combine denominators 

(1 - z)Di + ^[(1 - u)Di + u[xD2 + (1 - x)D^]] = p"^ - A„ (24) 
p" = p — k[\ — z{l — ux)] — q[\ — z{l — u{\ — x))], 
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A^^ zA- z{l - z)[k{l - ux) + q{l - u{l - x))]^ + (1 - z)m'^, 
^ ^ ^ ^ = 6 / dx du dzuz^- — ^— tt. 

DlD2^3^4 Jo Jo Jo (p"2-A^)4 

Note that A|^=i = A and pj'^^^ = p'. 

After the shift in the crossed diagram we have 

4/300 = -6 /' dx [' du [' dzuz' [ '^{\p"'C^^l, (25) 
Jo Jo Jo J I 6 

where 

^SoO ^ —9a09l30 + 5*0^ = --^ajoO' (^6) 

C'Soo = ^k'^gaogpoo^ - ^k^gaogpoa + '^k'^gaogoo - ^k'^gapa^ + ^k'^gapa (27) 

-^^5'a/3 + 'i-^{kq)gaogpoab - S{kq)gaogpoa - S{kq)gaogpob + ^{kq)gaogpo 
-12{kq)gapab + 6{kq)gapa + 6{kq)gapb - 3{kq)gap + Ak^kpa^ - Akak^a 
+2kakp - Sk^qogpoab + Ak^q^gpoa + 4A;«gofl'/30^ - 2/c„gofl'/3o - Sk^qog^oab 
+Ak(iqogaoa + Akf^qog^ob - 2A;^go5'aO + 8g^5'«o5'/30&^ - Sq'^gaogpob + ^.q^gaogpo 
-^q^gapb'^ + ^q^gapb - q^gap + ^qlgapb^ - ^qlgapb + ^qlg^p - Agaogpom^ 

+4:gapm'^: 

^apoo = '^k'^gaogpoa^ " ^k'^gaogpoO^ + '^k'^gaogpoo^ - k'^gapa^ + '^k'^gapa^ (28) 

-k^gapa'^ + 8k'^{kq)gaogpoa% - 4:k'^{kq)gaogpoa^ - \.'2k'^{kq)gao9po0^b 
+Qk'^{kq)gaogpoa^ + 4:k^{kq)gaogpoab - 2k'^{kq)gaogpoa - Ak'^{kq)gapa% 
+2k'^{kq)gapa^ + 6k'^{kq)gapaH - 3k'^{kq)gapa'^ - 2k'^{kq)gapab + k^{kq)gapa 
-4k'^kaqogpoa%+2k'^kaqogpoa^+^k'^kaqogpoa^b-3k'^kaqogpoa'^-2k'^kaqogpoO'b 
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+k'^kaqog(3oa - ik^k^qogaoa^b + 2k'^k/3qogaoa^ + Qk'^kpqogaoo^b - 3k'^kf^qogaoa^ 
-2k'^kpqogaoab + k'^k^qogaoa + Ak'^q^gaogpoa^'^ - Ak'^q^gaQgpoa^b 
+2k^q'^g^ogpoa'^ - Ak'^q^g^og^oab^ + Ak'^q'^gaogpoab - 2k'^q'^g^ogpoa 
+2k'^q^gaogpob'^ - 2k'^q^gaogpob - 2k'^q^gapa%'^ + 2k^q^gaf3a% - k'^q^gapo^ 
+2k^q'^gai3ab^ - 2k'^q^gapab + k^q^gapa - k^q'^gapb^ + k'^q^gapb + 2k'^qlgapa^ 
-'^k'^qlgapa - 4:k^gaogf}oa'^rn'^ + ^^k'^gaogpoam'^ - '^k'^gaogpom'^ + 2k'^g^i3a^m^ 
-2k'^gapam'^+k'^gapm'^+S{kqygaogpoa%'^-S{kqfgaog(3oa'^b-8{kqygaogpoab'^ 
+8{kqfgaogf3oab-A{kqfgapa^b^ + A{kqfgapa^b + A{kqfgaf3ab^-4:{kqfgapab 
-S{kq)kaqog/3oa^b^ + S{kq)kaqogpoa^b + 8{kq)kaqogpoab^ - S{kq)kaqogi3oab 
-9>{kq)kpqogaoa%^ + S{kq)kpqogaoa'^b + 8{kq)kpqogaoab'^ - S{kq)kpqogaoab 
+S{kq)q'^gaogpoab^ - 12{kq)q^gaogpoab'^ + 4:{kq)q'^gaogpoab - 4:{kq)q'^gaog/3ob^ 

+6{kq)q^gaogi3ob^ - 2{kq)q^gaQgpob - A{kq)q^gapab^ + 6{kq)q^gapab^ 
-2{kq)q'^gaf}ab+2{kq)q'^gaf)b^ - 'i{kq)q^gai3b'^ + {kq)q'^gaf}b-8{kq)gaogpoabm^ 
+4:{kq)gaogpoam'^ + A{kq)gaogpobm^ - 2{kq)gaogpom'^ + 4:{kq)gapabm'^ 
-2{kq)gai3am'^ - 2{kq)gagbm'^ + {kq)gafim^ + 2kakpq%'^ - 2kak/3q'^b 
+8kakpq^a%'^ — Skakpq^a^b — Skak/^q^ab^ + Skak/^q^ab — 2ka^kpm^ 
-Akaq^qogpoob^ + Qkaq^qogp^ab^ - 2kaq'^qogpoab+2kaq'^qogpob^ - ^kaq^qog^ob'^ 

+kaq^qogf3ob + Akaqogpoabm^ - 2kaqogi3Qam^ - 2kaqogf}obm^ + kaqogpom^ 
-Akpq^qogaoab^ + dk^q'^qogaoab'^ - 2k,3q'^qogaoab + 2k,3q'^qogaob^ - 3/c/3g^go5'ao&^ 
+kpq'^qogaob + Akpqogaoabm'^ - 2kpqogaoam'^ - 2kpqogaobm'^ + k^q^gaom^ 

+'^q'^gaogf}ob^ - ^^gaog/^ob^ + '^q^gcogpob'^ - q^gaf^b'^ + "^q^gapb^ - q^ga^b'^ 

-Aq^gaognob'^m'^ + Aq'^gaog^obm'^ - 2q'^gaogpom'^ + 2q^gapb'^m^ - 2q^gapbm'^ 

+q^gai3m^ - 2qlg^pm'^ + 2g^QgpQm'^ - gapm'^, 
a = 1 — z{l — ux), 
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b= 1 - z[l - u{l - x)]. 
Momentum integration leads to 

Ca/300 = (29) 

Again we use integration by parts over z to separate the ultraviolet di- 
vergent term 

dzz" log — = — log —\l 30 

A2 3 A2 

"1 dz 



L 



dz 



1 /•! 
-log — - / 
3 JO 



3A 



-{l-z-^){{^-m')-{l-2z)[k{l-ux) + q{l-u{l-x))Y]. 



Subtraction in the logarithmic term at zero momenta is performed by 
omitting the log as well as the constant term after integration by parts, as 
was done for the ladder diagram. Again this subtraction does not exhaust all 

necessary subtractions. The function C^^poo contains terms {—'igaogfso'm'^ + 

Agapfn^) and function Cj^^Q^ contains terms 2gaog/3om'^ — gap^'^ which also 
require subtractions. 

To perform subtractions this time it is sufficient to replace 

11 1 



A, A,(A;,g) A,(0,0) 
z{w? - A) + z{l - z)[k{l - ux) + q{l - u{l - x))f 



(31) 



m^A 



z 



Al ^ Al{k,q) A2(0,0) ^A,{k,q) A,(0,0)^^A,(A;,g) + A,(0,0)^ 
z{m'^ - A) + z(l - z)[k(l - ux) + q(l - u(l - x))]^ 



{m' + A,)- 



m 



Al 



z 



in the denominators of the terms just discussed. 
After subtraction we obtain 
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c. 



•sub 



(32) 



■a/300 ~ 



where 



;^{l)sub _ 
'.,-.11 — 



(33) 



+12A;^firaofl'/3o^^^^^(l - ^) + Sk'^gao9i3ovz^{ux) + 4:k'^gaogi3ov{ux) 
+4:k'^gaogf3oz'^ {uxf- SPgaogpoz"^ (ux) + 4:k^gao9poz^ - '^k'^9ao9i3oz^ {uxf{l-z) 
+Sk^9aogpoz^{ux){l- z) - ik^ gaog/soz^ {1 - z) -6k^ gaog^oz^ -4:k^ gao9i3oz{uxf 
+&k'^9ao9poz{ux)-Ak'^9ao9f}oZ + 4:k'^9ao9f}o{uxf{l -z)- 8k'^gao9/3o{ux){l - z) 

+4A;2^„o^/3o(l - ^) + l^k'^gapo^z^ - l^k'^gapaz^ - 12k'^gapvh^{l - z) 
—8k'^gapvz^{ux) — Ak^gaf3v{ux) — Ak'^gaf3z'^{uxY + 8k^gaf3z'^{ux) — Ak'^gafsz^ 
+Aegapz\uxf(l -z)- 8k:'g^Bz'^{ux){l - z) + A^gc^BZ^l - z) + Sk^g^pz^ 

+Ak^9af3z{uxf - 8k'^gapz{ux) + Ak'^gaf^z - Ak'^g^piuxfil - z) 
+8k'^ga(3{ux){l - z)- Ak'^gapil - z) - A8{kq)gaogi3oabz'^ + 24:{kq)gaogpoaz'^ 
+24:{kq)gaogi3obz'^ + 24:{kq)gaogpoSvz^{l - z) + 8{kq)gaogpoZ^{ux)[u{l - x)] 
-8{kq)gaogi3oz'^{ux) - 8{kq)gc,o9poz'^[u{l - x)] + 8{kq)gao9f3oZ^ 
-8{kq)gao9i3oz^{ux)[u{l - x)]{l - z) - 16{kq)gaogi3oz^{ux)[u{l - x)] 
+8{kq)gao9i3oz^{ux){l - z) + 8{kq)gaogf3oz^[u{l - x)]{l - z) 

-8(fcg)fi'aofi'/3o^^(l - z)- 12{kq)gao9fioz'^ - 8{kq)gaog^oziux)[u{l - x)] 
+8{kq)gaogf3oz{ux) + 8{kq)gaogf3oz[u{l - x)] - 8{kq)g 
+8{kq)gaogi3o{ux)[u{l - x)]{l - z) - 8{kq)gaogpo{ux)[u{l - x)] 
-8{kq)gaogpo{ux){l - z) - 8{kq)gaogpo[u{l - x)]{l - z) + 8{kq)gaogf3o{l - z) 
-\-2>Q{kq)gaf}ahz^ — 18{kq)gai3az^ — 18{kq)gaj3hz^ — 2A{kq)gai3Svz^{l — z) 
-8{kq)gapz'^{ux)[u{l - x)] + 8{kq)gaj3z'^{ux) + 8{kq)gai}z'^[u{l - x)] 
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-8{kq)gapz'^ + 8{kq)gai3Z'\ux)[u{l - x)]{l - z) + 16{kq)gapz^{ux)[u{l - x)] 
-8{kq)g^f)Z^{ux){l -z)- 8{kq)g^pz^[u{l - x)]{l - z) + 8{kq)g^pz^{l - z) 
+9{kq)ga/3Z^ + S{kq)gapz{ux) [u{l - x)] - 8{kq)gapz{ux) - 8{kq)gaf)z[u{l - x)] 

+8{kq)gapz - 8{kq)gai3{ux)[u{l - x)]{l - z) + 8{kq)gap{ux)[u{l - x)] 
+8{kq)gap{ux){l- z)^8{kq)gap[u{l-x)]{l- z)-8{kq)gaf3{l- z)-12kakpa^ z^ 
+12kak/3az'^ - Qkakpz'^ + 2Akaqogpoabz'^ - '^'^kaqogpoaz^ - 12kaqo9i3obz'^ 
+^kaqo9poZ^ + 2Ak0o9aoabz^ - 'i-'^kgqogaoaz^ - Uk^qogaobz^ + Bk^qogaoz^ 
-2Aq^gaQgpoh'^z^ + 2Aq'^gaogpobz^ + 12q'^gaogf3os'^z^{l - z) 
+8q^gaogpQSZ^[u{l - x)] + Aq^gaogiios[u{l - x)] + Aq^gaogpoz'^[u{l - x)f 
-8q^gaQgpoz'^[u{l - x)] + Aq^g^ogpoz'^ - Aq^gaogi}oZ^[u{l - x)f{l - z) 
+8q^gaQgpQZ^[u{l - x)]{l - z) - Aq'^gaogpoz^il - z) - Qq^gaog^QZ^ 
-Aq^gaQg(3Qz[u{l - x)f + 8q^gc,QgpQz[u{l - x)] - Aq^gaogpoZ 
+V5'ao5/3o['"(l - x)f{l- z) - 8q'^gao9f3o[u{l - x)]{l - z) + Aq^9ao9l3o{.'^ - ^) 
+I8q^g^pb''z'' - 18q^gapbz'' - 12q''g^f^sh^{l - z) - 8q''gapsz^[u{l - x)] 
-4q^gaps[u{l - x)] - Aq^gapz'^[u{l - x)f + 8q^gapz^[u{l - x)] - Aq^g^pz"^ 
+V^a/3^'[«(1 - x)Wl -z)- 8q^g^pz^[u{l - x)]{l - z) + Aq^g^^z^l - z) 
+^^9apz'^ + V5'a/3^[ii(l - 2;)]^ - 8q'^gafsz[u{l - x)] + Aq^g^pz 

- Vyc«/3 K 1 - ] ' ( 1 - -Z) + 8?'^a/3 K 1 - a:;) ] ( 1 - -2) - V ( 1 - -Z) - 1 2?o ^a/3&' 

+ ^'^<ll9alibz'^ - Qql9af3Z^: 

^alot = '^k^9ao9i3oa'^ - 4^k'^9ao9poa^ + '^k^9ao9poa^ - 2A;^5'ao5'/3o^^^^^(l - z)^ 

(34) 

-Ak'^gaog(3ov^z'^{ux){l-z)-2k^gaogi3ov'^z'^{uxf-k'^gai3a'^+2k'^gai3a^-k^gai3a^ 
+k^g^pv^z^{l - zf + 2k^g^pv^z''{ux){l - z) + k^g^pv''z'{uxf 
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+8k'^{kq)gaog^Qa% - Ak'^{kq)gaogpoa? - 12k^{kq)gaQgpoa?h + &k'^{kq)gaogpoo? 
+Ak'^{kq)gaogf)oab - 2k'^{kq)gaogf}oa - 8k'^{kq)gaogf)oSv^z'^{l - zf 
-%k^{kq)gaQgi3QSV^z^{ux){l - z) + 8k'^{kq)gaogpov'^z'^{ux)[u{l - x)]{l - z) 
+8k'^{kq)gaogpovz'^{uxf[u{l - x)] - Ak^{kq)gapa% + 2k^{kq)gapa^ 
+6k'^{kq)gaf3a^b — 3k'^{kq)ga/3a^ — 2k'^{kq)gai3ab + k'^{kq)gapa 
+4:k'^{kq)gai3sv^z'^{l - z)'^ + Ak'^{kq)gapsv'^z^{ux){l - z) 
-Ak^{kq)gai3V^z^{ux)[u{l - x)](l - z) - Ak'^{kq)gai3vz^{uxf[u{l - x)\ 
-4:Pkaqogi3oa^b+2Pk„qogpoa'^ + 6k^k^qogpoa^b-3k^kaqogi3oa^-2k^kaqogpoab 
+k'^kaqogi3oa - Ak^kpqogaoa^b + 2k^kpqogaoa^ + Gk'^kpqogaoo^b - Sk'^kgqogaoo^ 
-2k^kf)qogaoab + k^k^qogaoa + Ak'^q^gaogj3oa%'^ - ^^k'^q'^Qaogf^oa'^b 
+2k'^q^g^ogpoa^ - Ak'^q^g^ogpoab^ + Ak^q'^gaogpoab - 2k^q^g^ogpoa 

+2k^q^go,oggob^ - 2k^q^gaoggob - ik^q^gaog^os^v^z^il - zf 
-Ak'^q^gaQgpQS^vz^{ux){l - z) - Ak^q^gaogposv^z^[u{l - x)]{l - z) 
-Ak^q^gaQgi3QSVz^{ux)[u{l - x)] - 2k^q^gai3a%'^ + 2k^q^gai3a% - k^q^gapo^ 
+2k^q^g„pab^ - 2k^q'^g^pab + k'^q^g^pa - k'^q^g^nb^ + k^q^gaf^b 
+2k^q^g^ps^v^z'^{l - zf + 2k^q^g^gs^vz\ux){l - z) 
+2k'^ q'^ ga/3sv'^ z'^[u{l — x)]{l — z) + 2k'^q'^gai3svz'^{ux)[u{l — x)] + 2k'^q^gai3C? 
-2k'^qlgai}a - Ak'^gaognoa^m'^ + Ak'^gaogpoaw? + 2k'^gaogi3orn^v'^z'^{l - z) 
+2k'^gaogi3om\'^z{l - zf + 2k'^gaQgf)Qm^v^z{l - z) ^ 2k'^gaogi3om'^vz'^{ux) 
+2k'^gaog8om'^vz{ux){l- z) ^-2k'^ g„QgfiQm^vz{ux)-2k^ g^^gf^^m^ ^2k^ g^^a^vr? 
-2k^gapam^ - k^ g^firri^v^z^iX -z)- k^ ga,firr?v^z{\ - zf - k^ga.firri}v^z{\ - z) 
—k'^gap'>TT''^'vz'^{ux) — k'^ gai3rn^vz{ux){l — z) — k'^gai3nn?vz{ux) + k'^gapm^ 
+8{kqfgaogf3oaH^ - 8{kqfgaogpoa% - 8{kqfgaogf3Qab'^ + 8{kqfgaogpoab 
-8{kqfgao9pos'^v'^z'^{l - zf + 16{kqfgaogi3oSvz'^{ux)[u{l - x)]{l - z) 
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-8{kqygaogi3oz^{uxY[u{l-x)Y-4{kqfga^a^b^+4:{kqYgo,i3a^b+4{kqYgai3ab'^ 
-A{kqfg^pah^A{kqfg^ps'^v'^z^{l - zf - 8{kqfgc,i3svz^{ux)[u{l - x)]{l - z) 
+A{kqfgapz^{uxf[u{l - x)f - 8{kq)kaqogpoa%'^ + 9>{kq)kaqogpQa% 
+8{kq)kaqogpoab'^ - 8{kq)kaqog poab - 8{kq)kpqogaoa^b'^ + 8{kq)kpqQgaQa% 
+8{kq)kpqogMab'^ - 8{kq)kfsqogaoab + 8{kq)q'^gaogf3oab^ - l2{kq)q^gaogj3oab'^ 
+4:{kq)q'^gaogfSoab - 4:{kq)q'^gaogpob^ + Q{kq)q'^gaogpob'^ - 2(/cg)g^5'ao5'/3o& 
-8{kq)q'^gaog(3oS^vz'^{l - zf - 8{kq)q'^gaogpos'^vz'^[u{l - x)]{l - z) 
+8{kq)q^gaogi3oS^z^{ux)[u{l - x)]{l - z) + 8{kq)q^gaogpoSz'^{ux)[u{l - x)f 
-A{kq)q^g^pab^ + &{kq)q^gapab'^ - 2{kq)q^ gapab + 2{kq)q^ g^gb^ - 2>{kq)q^ gajib'^ 
H^q^Qapb + A{kq)q^gaps^vz'^{l - zf + A{kq)q^g^ps^vz'^[u{l - x)](l - z) 
-A{kq)q^gai}S^z'^{ux)[u{l - x)\{l - z) - A{kq)q^gai3sz^{ux)[u{l - x)f 
-8{kq)gaogi3oabm^ + A{kq)gaQgpQam^ + A{kq)gaogpQbm'^ 
+A{kq)gaogpom'^svz^{l - z) + A{kq)gaog()om^svz{l - zf 
+A{kq)gaogpom'^svz{l - z) - A{kq)gaogpom'^z'^{ux)[u{l - x)] 
-^kq)gaogi3om^z{ux)[u{l - x)]{l - z) - A{kq)gaogpom^z{ux)[u{l - x)] 
-2{kq)gaogpom'^ + A{kq)g,^0abm^ - 2{kq)ga(3am'^ - 2{kq)g,,3bm^ 
—2{kq)gai3nn?'svz^{l — z) — 2{kq)ga(3m^svz{l — zf — 2{kq)gai3m^svz{l — z) 

■^2{kq)gai3rn^ z^ {ux)\u{l — x)] + 2{kq)gai3nn? z{ux)\a(l — x)](l — z) 
+2{kq)gaprn'^z{ux)[u{l—x)] + {kq)gai3rn'^+2kakpq^b'^—2kaki3q'^b+8kakpqQa'^b'^ 

—8kak/3qla^b — 8kakpq^ab'^ + 8kakgq^ab — 2kak/sm'^ — Ak^q^q^gp^ab^ 
+^K(f%gi3Qab'^ - '^Kiftlogpoab + 2A;„g^go5'^o&^ - '^Kif^gpob'^ + Kq'^qo9/3ob 
+Akaqogf3oabm'^ - 2kaqog/3oam'^ - ^kaqog/sobm^ + kaqogpom^ - Akpq^qogaoab^ 
+6kpq'^qogaoab'^ - 2kpq'^qogaoab + 2kpq^qogaob^ - ^kpq'^qogaob'^ + k^q^q^gaob 
+Akpqogaoabm'^ - 2A;^Qo5'aoa^^ - '^k^ogaobm'^ + k^qogaom'^ + '^q'^gaogpob'^ 
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-M'^9ao9(5oS^z^ [u{l-x)]{l-z)- g^ogpQS^ [u{l-x)f~ q^gaph'^ + ^q^Qaph^ 
-q^gaph^+q^gaps^z^{.l-zf+2q^g^ps^z'^[u{l-x)]{l-z)+q^g^ps^z^[u{l-x)f 

-4^q^gaog(3ob^rn^ + 4:q^gaogpobm'^ + 2q^gaogpom^s^z'^{l - z) 
+2q^gaQgpom^s^z{l - zf + 2q'^gaogpom^s'^z{l - z) + 2q'^gaogf3om^sz'^[u{l - x)] 
+2q^gaQgpom^sz[u{l - x)](l - z) + 2q^gaogpom^sz[u{l - x)] - 2q^gaogpom^ 

, r, 2 1.2 2 o2 7,2 2 222/1 \ 2 22/1 \2 

+2q gapo m -2q gapom - q ga^m s z [I - z) - q gapm s z[l - z) 
-q^gapm^s^z{l - z) - q^gapm^sz^[u{l - x)] - q^gapm^sz[u{l - x)]{l - z) 
-q^gapm^sz[u{l - x)] + q^gapm^ - 2qlgapm'^. 

3.3 LOW FREQUENCY LIMIT 

We have obtained in eq.(|20|) and eq.(p^) explicit expressions for the con- 
tributions to the hght by hght scattering tensor induced by the ladder and 
crossed diagrams, respectively. These expressions have the form of the two- 
and three-dimensional integrals over the Feynman parameters. A very impor- 
tant check of the correctness of all preceding calculations may be performed 
now. Namely, we are going to compare the low-frequency limit of the ex- 
pressions above for the light by light scattering tensor with the well known 
expression given by the Euler-Heisenberg Lagrangian. Leading terms in the 
expansion of the tensor in eq. over k and q should coincide with the Euler- 
Heisenberg Lagrangian. Naively, it is far from evident that the expansion of 
the integrals over Feynman parameters obtained above will start with terms 
at least quartic in momenta and that these leading terms will be at least 
quadratic simultaneously in q and k. 

Performing the low frequency expansion in eq.(pOD and eq.(^) we obtain 
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1 7 4 1 

1 13 38 4 

1 3 7 1 

+ 12^^^^'^°^"° + 30^^^^°^"°"^' + lij^"^'^"^^" + 30^^"^°^^°^'- 

Note that some undesirable terms containing too small a number of fac- 
tors of or g emerged in the expansion. One should not be too disappointed 
at this stage since we have performed above a transformation of the light 
by light scattering tensor taking into account that further integration over k 
kills all odd factors in k. It is easy to see that all unwanted terms are odd in 
k and, hence, they disappear after the integration over k. 

Really, it is not difficult to calculate second ladder diagram for the light 
by light scattering separately and to obtain a low frequency asymptote which 
does not contain unwanted terms at all 



45m4 ^ ' 9m4 " ""'^ A5m^ 

1 7 , , o 1 



,2 



-ikq)kaqo9(30 - T^k^k^q + T^^kpk^q^ 



It is easy to see that this last expression coincides exactly with the one 
corresponding to the Euler-Heisenberg Lagrangian (for the Euler-Heisenberg 
expression see, e.g. [Q). 

4 CALCULATION OF THE CONTRIBUTION 
TO THE LAMB SHIFT INDUCED BY 
THE LIGHT BY LIGHT SCATTERING 
INSERTION 
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4.1 LADDER DIAGRAM 

We start the calculation of the contribution to the energy shift induced by 
the ladder diagrams in Fig. 1 by multiplication of the respective part of the 
light by light scattering tensor by the numerator of the electron factor in 
eq.(0)0. Thus we obtain 

Li = 2Aai3C^am (37) 

where 

= —2k'^qoUvx{l — u) — Qk'^qoUv{ux) + Qk'^ qou{ux)'^ — Ak'^qou{ux) (38) 

—k'^qoVx{l — u) + k'^qox{l — u){ux) + 2A;^go(l — u){uxY — Ak'^mu{uxY 
+2k'^mu{ux) + ?)2{kq)qQu{ux)[u{l — x)] — 2{kq)qQu{ux) — 2{kq)qQu[u{l — x)] 
—2{kq)qou + 4:{kq)qoxl{l — u){ux) + 8(/cg)go(l — u){ux)[u{l — x)] 
—^{kq)mu{ux)[u{l — x)] + 2{kq)mu{ux) + 2{kq)mu[u{l — x)] + 2{kq)mu 
—Qq^ qQSu\u{l—x)]— q^ qQSx\{\—u) —2q^ qQs{\—u) \u(l—x)] + \Qq^ qQu\u{l—xyf' 
-2q^qQu[u{l - x)] + q^q^xlil - u){u{\ - x)\ + 2g^go(l - u){u(\ - x)^ 
-A(^mu\u{\ - x)f + 2(^mu\u{\ - x)\ - %(^qu\u{\ - x)]^ - 4goM[M(l - x)\ 
\%(^mu\u{\ — x)f + 4gQmM[n(l — x)], 

= k^qQv'^{uxY — k^mv'^{uxY + k'^m{uxY — k^m(uxY (39) 

—4:k'^{kq)qov{uxy[u{l — x)] — 4A;^(A;g)go(wx)^['u(l — x)] + k'^{kq)qQ{ux)^ 

+3k'^{kq)qo(uxY[u{l — x)] — k'^{kq)qo(uxY + Ak'^ {kq)mv{uxY[u{l — x)] 

+Ak'^ (kq)m(ux)'^[u{l — x)] — k'^ (kq)m(ux)^ — 3k'^ {kq)m(uxY[u{l — x)] 

+k'^ {kq)m{uxY + 2k'^q^qQSv{ux)\u{\ — x)] — k'^q'^qo{uxY[u{l — x)]^ 

+k^q'^qo{uxy[u{l - x)] + 2k^q^qo{ux)[u{l - x)]^ - 2k'^q'^qo{ux)[u{l - x)] 
■^Additional factor 2 is inserted to take into account both ladder diagrams in Fig.l 
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—2k'^q'^msv{ux)[u{l — x)] + 2k'^q'^m{ux)'^[u{l — x)]^ — /c^g^m(-ux)^['u(l — a:)] 

-k^q'^m{ux)[u{l - x)]'^ + 2k^q'^m{ux)[u{l - x)] - 2k'^qQ{ux)'^[u{l - x)]'^ 

+2k'^q^{uxf[u{l - x)] - 2k'^qlm{uxf[u{l - x)f - 2k'^qQm{uxy[u{l - x)] 

—2k'^qom'^v{ux) + k'^qom^{uxY — 2k'^ qom^ (ux) + 2k'^m^v{ux) — 2k'^m^{uxy 

+k^m^{ux) - 4:{kqfqo{uxf[u{l - x)]^ + 2{kqfqo{uxf[u{l - x)] 

+2{kqfqo{ux)[u{l - x)]^ - 2{kqfm{uxf[u{l - x)] - 2{kqfm{ux)[u{l - x)]^ 

-4:{kq)q'^qos{ux)[u{l - x)]'^ - 6{kq)q'^qo{ux)[u{l - x)]^ 

+S{kq)q''qo{ux)[u{l - x)f + {kq)q\oHl - x)]^ - {kq)q''qo[u{l - x)]' 

+4:{kq)q'^ms{ux)[u{l — x)]'^ + 4:{kq)q'^m{ux)[u{l — x)]^ 

—3{kq)q^m{ux)[u{l — x)]'^ — {kq)q^m[u{l — x)]^ + {kq)q^m[u{l — x)]'^ 

+A{kq)qo{ux)[u{l-x)]^-A(kq)qlm{ux)[u(l-x)]^ + 10{kq)qom'^(ux)[u{l-x)] 

— {kq)qo'm'^{ux) — {kq)qo'm'^[u{l — x)] — {kq)qQW? — '&{kq)m^ {ux)[u{l — x)] 

+ {kq)m^{ux) + {kq)m^[u{l ~x)] + {kq)m^ + q'^qos'^[u{l — x)]"^ — q'^qo[u{l — x)]'^ 

+q'^qo[u{l — x)]^ — g^ms^[-u(l — x)]^ + g*^m[-u(l — x)]'^ — q'^m[u{l — x)]^ 

+2q'^ql[u{l - x)f - 2q'^q^[u{l - x)]^ - 2q'^q^m[u{l - x)]^ + 2q'^qlm[u{l - x)]^ 

-2q^qom^s[u{l-x)] +2q^qom'^[u(l-x)]'^-q^qom'^[u{l-x)] +2q^m^s[u(l-x)] 

-2gW[M(l - x)f + - x)] - 2qlm'^[u{l - x)f - 2q^m^[u{l - x)] 

+2qlm^[u{l - x)]'^ + 2qlm^[u{l - x)]. 

Next we are going to calculate the upper loop momentum integral in 
Fig. la, 6. Due to two photon propagators this integral contains a factor q'^ 
in the denominator and one may worry about its convergence at small loop 
momenta. We have checked in the previous section that the light by light 
scattering tensor supplies at least two powers of loop momentum q in the low- 
frequency region making the loop integration completely safe. However, the 
integral over the upper loop momentum q will be taken prior to integration 
over the Feynman parameters in the light by light scattering loop and in 
this case the problem of poor low q behavior of the integrand resurrects. 
In order to make all intermediate integrations infrared safe we temporarily 
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introduce an intermediate photon mass A for the photons in the upper loop. 
All calculations both for the ladder and for the crossed diagrams will be 
performed with this nonvanishing photon mass. The final expression for 
the energy shift will admit the limit of vanishing photon mass and will be 
convergent even in this limit. 

Let us combine both electron denominators and photon denominators in 
the upper loop in Fig.l. To facilitate transformations we introduce new a 
notation 

--{q'^ + ak'^ -2p{kq) --fm^), (40) 
7 



where 



x(l — ux) 

a — 



{I - x)[i - u{i - x)y 

ux 



i-u{i-xy 
1 



7 = 



u{l - x)[l - u{l - x)]' 



Then 

(1 - t)[{l - y){q' - X') + y{q' + 2mgo)] + t{-^A) = q'' - dl (41) 
where 

q' = q + y{l — t)m — j3tk, 
dl = -kH{a - (3H) + m''[y\l - tf + 7^] + \\l -t){l-y). 



Hence, 

1 



-,,j\y(l-y)j\t(l-tT--^^, (42) 



(g2 - A2)2(g2 + 2mgo)A 'Jo "'Jo^ ' {q'^ - dl) 



1 



24f/*(l-,)/d«(l-i)^ 1 



(g2 - A2)2(g2 + 2m?o)A2 ' Jo ^'Jo ' ' {(l''^-dlY 
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Next euclidean rotation is performed followed by calculation of the inte- 
gral over q'^ for the contribution of the ladder diagram to the energy shift 
in eq.(|lOD. 

m = 3 m{-^f—a\Zafm{-^f dx (43) 



r du-f 1^ dy{l-y) dt{l - tf dx du f dq'^q'^ 
Jo Jo Jo Jo Jo J 

F[{q = q' - y{l - i)m + ptk) 



{{l-u) 



(g'2 + dlY 



+87^1 - u)t i^i^^^ >• 

Explicit expressions for the numerators in the integrand (after Wick ro- 
tation) have the form 

2F[{q') = -q''^mP[ + 2PmP^ + 2m^[y{l - t)fPl (44) 

where 

Pi = 18su[u{l - x)][y{l - t)] + 3s{l - x){l - u)[y{l - t)] (45) 
+6s(l - u)[u{l - x)][y{l -t)] - 18u[u{l - x)]'^[y{l - t)] - 4u[u{l - x)]'^ 
+12u[u{l - x)][y{l - t)] + 6m[m(1 - x)] - 3(1 - x)(l - u)[u{l - x)][y{l - t)] 

-6{l-uMl-x)]'[y{l-t)], 



P^ = 6su[uil - x)Ml - mptf + 3(1 - x)(l - uMl - t)]iPtY (46) 

+2s(l - u) [u{l - x)] [y{l - t)]{(3tf + 2uvx{l - u) [y{l - 1)] + 6uv{ux) [y{l - 1)] 

-6m(mx)^[?/(1 - t)] - Au{ux)^ - 32u{ux)[u{l - x)][y{l - t)]{(3t) 
-8u{ux) [u{l- x)]{l3t) + 2u{ux) [y{l- t)]{l3t) + Au{ux) [y{l-t)]+ 2u{ux) (fSt) 
+2u{ux) - 10u[u{l - x)]^[y(l - t)]{ptf - 4:u[u{l - x)Y{ptf 

+2u[u{l - x)][y{l - t)]{ptf + 2u[u{l - x)][y{l - t)]{pt) + 2u[u{l - x)]{ptf 
^Subscript E is omitted below. 
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+2u[u{l - x)]{pt) + 2u[y{l - t)]{pt) + 2u{(3t) + vx{l - u)[y{l - t)] 

-x{l - u){ux)[y{l - t)] - 4(1 - a;)(l - u){ux)[y{l - t)]{pt) 
_(1 - x)(l - u)[u{l - x)Ml - t)](/3t)2 - 2(1 - u){uxy[y{l - t)] 
-8(1 - u)iux)[u{l - x)Ml - t)]m - 2(1 - uMl - x)ny{l - t)]{(3tf, 

Pi = Qsu[u{l - x)] [y{l -t)]+s{l-x){l-u) [y{l - t)] (47) 
+2s(l - n)[n(l - x)][y{l - t)] - 2u[u{l - x)]^[yil - t)] + 4m[u(1 - x)f 
+6u[u{l - x)][y{l - t)] + 6u[u{l - x)] - (1 - x)(l - u)[u{l - x)][y{l - t)] 

-2{l-u)[u{l-x)ny{l-t)]. 
The second term in the integrand in eq.(^) has the following form 

2F^ = ^q'^m[u{l - x)]'^T^ - q'^k'^m[u{l - x)]T/ - g'W[M(l - x)]T^ (48) 
+2k^mT^ + 2k'^m^Tl + 2m^[M(l - x)][y{l - t)fTl 

where 

Tl = -12s2[y(l - t)] - 6s' + 2[u{l - x)]'[y{l - t)] (49) 
+3[u(l - x)]' - 2[u{l - x)][y{l - t)] - 3[u(l - x)], 

Tl = -8s'[u{l - x)][y{l - t)]{f3tf - 6s'[u{l - x)]{(3tf (50) 
-Qsv{ux)[y{l - t)] - Asv{ux) + lQs{ux)[u{l - x)][y{l - t)]{l3t) 
+12s(mx)[m(1 - x)]{pt) + 8{uxY[u{l - x)][y{l - t)] + 3{uxf[u{l - x)] 
-7{uxf[y{l - t)] - 4:{uxf + 18(na;)[n(l - x)f[y{l - t)]{pt) 

+ 10(Mx)[M(l-x)]2(/3t)-12(Mx)[M(l-x)][?/(l-t)](/3t)-7(Mx)[M(l-x)][y(l-t)] 

-9{ux)[u{l - x)](/3t) - 3(mx)[m(1 - x)\ + 6(mx)[?/(1 - t)\ + 4(mx) 
+5[«(1 - x)f\y{\ - t)\{{itf + 5[n(l - x)\\{itf - 5[n(l - x)f\y{\ - t)\{(itf 
-4[m(1 - x)f{y(\ - t)](/3t) - 5[m(1 - xf{^tf - 3[m(1 - x)\'{i5t) 
+4[m(1 - x)\{y(\ - t)](/3t) + 3[n(l - x)](/3t), 
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= -8s'[u{l - x)] [yil - t)f - Qs'[uil - x)] [yil - t)f (51) 
+6s[y{l - t)] + 4s - 5[m(1 - x)f[y{l - t)f - 3[u{l - x)f[y{l - t)]^ 
+5[u{l - x)ny{l - t)]' + 3[u{l - x)ny{l - t)]' - 3[u{l - - t)] 

-3[u(l-x)]+6[y(l-t)]+3, 

n = -^'[^^(1 - x)ny{i - t)]{ptr - s'[u{i - x)npty (52) 

-2sv{ux)[u{l - x)][y{l - t)]{ptf - 2sv{ux)[u{l - x)]{ptf 

+4s{ux)[u{l - x)]'^[y{l - t)]{(3tf + 4s{ux)[u{l - x)]\(3t)^ 
-v^{uxf[y{l - t)] - v^{uxf + 4v{uxf[u{l - x)][y{l - t)]{l3t) 
+4:v{uxf[u{l - x)]{pt) + {ux)^ + 4:{uxf[u{l - x)][y{l - t)]{pt) 
+A(uxf[u(l - x)](pt) - {uxf[y(l - t)]((3t) - (ux)^{pt) - (ux)^ 

+5{uxf[u{l - x)]^[y{l - t)]{(3tf + 2{uxf[u{l - x)]\(3t)^ 
-3{uxf[u{l - x)Ml - t)]iPtf - 3iuxy[uil - x)Ml - t)]{(3t) 
-3{uxf[u{l - x)]{(3tf - 3{uxf[u{l - x)]{(3t) + {uxf[y{l - t)]{(3t) 
+(uxy((3t) + 6(ux)[u(l - x)]^[y(l - t)]((3t)' + A(ux)[u(l - x)]'((3t)' 
-3{ux)[u{l - x)ny{l - t)mf 4(u:r)[n(l - x)]'^! - 
-3{ux) [u{l- x)]\(3tf -3{ux) [u{l - x)]\(3t)^ + 2{ux) [u{l- x)] [y{l - t)]{(3tf 
+2{ux)[u{l - x)]{ptf + [u{l - xmi - t)]{pt)^ + [u{l - x)]\pt)^ 

-[u{i - x)ny{i - t)]m' - [u(i - x)ny{i - t)mr - hi - x)]w 

-[a{l - x)]%ptf + HI - x)nyil - t)]{ptf + [«(1 - x)]'{ptf, 

Ti = -2sM^ - ^)fW - t)?m' - 2«'[«(1 - - t)fm' (53) 

-2sv{ux)H'^ - x)][y{l - t)f - 2sv{ux)H'^ - x)][y{l - t)f 
+As{ux)H^ - x)f[y{l - t)f{(5t) + As{ux)H^ - x)f[y{l - t)f{(5t) 
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+2s[u{l - x)][y{l - t)]{ptY + 2s[m(1 - x)]{ptf + 2v{ux)[y{l - t)] + 2v{ux) 
+3{uxy[u{l - x)ny{l - t)]' - 3{uxnu{l - x)][y{l - t)]' 
-3iux)'[u{l - x)][yil - t)]' - iux)'[yil - t)] - 2iux)' 
+2{ux)[u{l - x)f[y{l - t)fm - 3iux)[u{l - x)]'[yil - t)fm 
-2{uxMl - x)ny{l - t)f - 3{uxMl - x)ny{l - t)]\pt) 
-MKl-x)]^[y(l-i)]2+2MKl-a;)][y(l-i)]3+2MKl-a;)][y(l-i)]^ 
-10{ux)[u{l - x)][y{l - t)]{pt) - 8{ux)[u{l - x)]{pt) + {ux)[y{l - t)]{pt) 
+2{ux)[y{l - t)] + {ux){(3t) + (ux) - [u{l - x)f[y{l - t)]^(/3t) 
-[«(l-a;)]3b(l-t)]2(/3t) + [«(l-x)]2b(l-t)]='(/5t) + [«(l-x)]2[l/(l-t)]2(/3t) 

-2[u{i - x)ny{i - mptf - 2[u{i - x)nptr + [u{i - x)][y{i - mptf 

n = -s'[u{l - x)][y{l - t)f - s'[u{l - x)][y{l - t)]' (54) 
+2s[y(l - t)] + 2s- [u(l - x)ny(l - t)]' - [u(l - x^l - t)]' 
+[u{l - x)ny{l - t)f + [uil - x)nyil - t)]' + 3[|/(1 - t)] + 3. 



After integration we obtain 




(55) 



mPl m^[y{l-t)]'^P^ 4 m[u{l - x)]'^uTl 
"3d|" ^ 3d| ^ 3^^ 2d| 




""i , m^\u(l — x)]u\t 




}■ 
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Integration over |k| may easily be performed since it is essentially one- 
dimensional due to spacelike nature of vector k. 

It is convenient to introduce the following notation 

dl = kH{a - f3H) + m2[y2(i _ tf + 7t] + - - y) (56) 



where 



a;,, = 



p = t{a-l3H), (57) 
m2[y2(i _ tf + 7t] + A2(l - t)(l - y) 



Integration over |k| in eq.(^) is impeded by the apparent infrared di- 
vergence, which is connected with the apparent constant asymptote of the 
electron factor. However, one expects that the leading low frequency term 
in this asymptote is proportional to k^ (if this is not the case the graphs 
under consideration would produce contributions of even lower order in Za 
which is well known to be the wrong conclusion). Surely, only the total 
electron-line factor should be proportional k^ in the low frequency limit and 
not the part taken into account in eq. (|55|) . These considerations give another 
opportunity to check the validity of all transformations above. One has to 
calculate all parts of the electron-line factor and then to calculate its low 
frequency asymptote. If it turns out that the electron-line factor vanishes 
as k^ when k is small then there are really no difficulties in the treatment 
of integrals of the sort contained in eq.(^); one has simply to perform sub- 
traction of the leading low frequency terms in the electron-line factor since 
these terms cancel in any case in the total electron-line factor. We will check 
the disappearance of the apparent leading term in the asymptote of the total 
electron factor in the next Section. 

Subtraction of the low frequency part is necessary in the last five terms 

of the integrand in eq.(^). The numerators of these terms are independent 

^We call here electron factor the product of the electron factor defined above and the 
light by light scattering tensor. Normalization of this electron factor differs from the one 
used in our previous works ^ on the contributions to the Lamb shift; in terms of the 
old normalization the asymptote we are discussing now is proportional to k^. 
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of momentum k and subtraction may be easily performed with the help of 
the identities 

1 1 k2 k2 



11 k2 k2 



(58) 



4(k) 4(0) pcoidlik) coldlik) 

k2 k^ 



- p^utik^ + ui) p^uiik^ + uiy^ 

11 k2 k2 k2 



4(k) 4(0) p'uldlik) putdUk) .;f4(k) 

k2 k2 k^ 



p3cu6(k2 + Lul) pM(k2 + ^!)' P'^!(k2 + cuf )3 • 
Substituting these expressions in eq.(^) we obtain 

A^r^ = "'(^f m(I^)3lg f'dx (59) 

Tm-^ m Jo 

rl rl rl roo 

/ du{l-u)-f dy{l~y) dt{l - tf d\k\ 
Jo Jo Jo Jo 



47* r ^,,2..rri ^3..rri m'Hl - x)]u[y{l - t)fTl 



^ 3p3(k2 + ^^)3 [-^"A^^2-^"^3 



[47t - mP^ + 2-ftm[u{l - x)]uT[ - 



1 M,,.^«^2 , o,.^L,.i .M.,^/ m3[y(l-t)]2pj 



3p2(k2 + cuf)2^ ' p 2 ' I I y n 1 

^^^^^ m3[n(l-x)]nT^ _ ^^^ nV>[u{l - x)]u[y{l - t)fTi 



3pa;2(k2 + c.2)^ 1 pc.2 

XKl-3:)]^r^ , , m^Kl-x)Hl/(l-t)]2T^ ^^ 
pui p^ul 

Integrating next over |k| we obtain 
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A^r^ = "'(^f f'dx (60) 

vrn-^ m n Jo 

1 rl rl 

du{l-u)-f dy{l-y) dt{l - tf 
Jo Jo 

7^ r ^,,2„,^/ ^3„,^/ m5[«(l-x)]iz[y(l-t)]2T' 



^Ap^i^l^ A 2 3 ^ 



A 



1 muTi , o . r .1 M ^/ m3[y(l-t)]2p^ 



A 

^^^^m^Kl-x)]^ _ ^^^ m^[u{l-x)]u[y{l-t)fTi 

OP^A /^A 

P^A P^^A 

4.2 CROSSED DIAGRAM 

Consideration of the crossed diagram contribution to the Lamb shift follows 
the same lines as in the case of the ladder diagram. First we multiply the 
respective part of the light by light scattering tensor by the numerator of the 
electron factor in eq.(^ and obtain 

Lc = AafsC^foQ (61) 



where 

= -Sk^qoa^z^ + SPqoaz^ + 6k'^qoV^z^{l - z) (62) 
^^k^qoVZ^{ux) + 2k^qov{ux) + 2k^qoZ^{uxf - Ak^qQZ*{ux) + 2k^qoz'^ 
-2k'^qoZ^{l - z){ux)^ + Ak^qoZ^{l - z){ux) - 2k^qoZ^{l - z) - k^q^z^ 
-2k^qQz{uxf + Ak^q^ziux) - 2^^50-2 + 2A;^go(l - z){ux)^ - 4A;^go(l - z){ux) 
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+2A;^go(l - z) - 2k^a^mz^ + 2Pamz^ - k^mz^ - 12{kq)qoabz^ + 6{kq)qoaz'^ 
+6{kq)qobz^ + 12{kq)qosvz^{l — z) -\- A{kq)qQz'^{ux)\u(l — x)\ — A{kq)qQz'^{ux) 
-A{kq)qoz'^[u{l - x)] + A{kq)qoZ^ - A{kq)qoZ^{l - z){ux)[u{l - x)] 
+A{kq)qQZ''{l - z){ux) + 4(A;g)go-2^(l - z)[u{l - x)] - A{kq)qQZ^{l - z) 
—S{kq)qQZ^ {ux)[u{\ — x)] — 3{kq)qoz'^ — 4:{kq)qoz{ux)[u{l —x)] +4:{kq)qoz{ux) 
+4:{kq)qoz[u{l — x)] — 4{kq)qoZ + 4{kq)qo{l — z){ux)[u{l — x)] 
-A(kq)qo(l - z)(ux) - A(kq)qo(l - z)[u(l - x)] + A(kq)qo(l - z) 
—A{kq)qQ{ux)[u{l — x)\ — A{kq)ahmz^ + 2{kq)amz^ + 2{kq)hmz^ — {kq)mz^ 
-lOq^qob'^z'^ + lOq^qobz^ + Gq^qos'^z^l - z) + Aq'^qoSZ^u{l - x)] 
+2q\os[uil - x)] + 2q^qoZ^[uil - x)]^ - Aq^qoZ^[u{l - x)] + 2q^qoZ^ 
-2q^qoz\l - z) [u{l - x)]^ + 4q\z\l - z) [u{l - x)] - 2q^qoz\l -z)- 2q^qoZ^ 

-2q''qoz[u{l - x)f + ^^ozHl - x)] - 2q^qQZ + 252^0(1 - z)[u{l - x)f 
-4g^go(l -z)[u{l-x)] + 2q^qQ{l-z)- 2q'^b^mz^ + 2q%mz'^ - q^mz^ + Aq^b'^z'^ 
-4.qlbz^ + 2qlz^ - iq^b'^mz'^ + Aqlbmz^ - 2qlmz^ , 

= 2k%a^ - Ak^qoa^ + 2k^qoa^ - 2k^qoV^z\l - zf (63) 
—4:k^qQV^z'^{l — z){ux) — 2k'^qQV^z^{uxY + k'^a'^m — 2k'^a?m + k^a^m 
-k'^mv'^z^il - zY - 2k^mv^z'^{l - z){ux) - k'^mv'^ z^ {uxY + 4:k'^{kq)qoa^b 
—2k^{kq)qoa^ — 6k'^{kq)qoa'^b + 3k^{kq)qoa'^ + 2k'^{kq)qoab — k'^{kq)qoa 
-8k'^{kq)qoSV^z'^{l - zf - 8k^{kq)qoSv'^z^{l - z){ux) 
+8k^{kq)qoV^z^{l - z){ux)[u{l - x)] + 8k^{kq)qoVz'^{uxf[u{l - x)] 
+4:k'^{kq)a%m - 2k'^{kq)a^m - 6k'^{kq)a'^bm + 3k'^{kq)a'^m + 2k'^{kq)abm 
-k^{kq)am - Ak^ {kq)msv^ z"^ {1 - zf - Ak'^{kq)msv^z'^{l - z){ux) 
+4fc^(/cg)mw^z^(l — z){ux)[u{l — x)] + Ak'^{kq)mvz'^{uxY[u{l — x)] 
+4:k'^q^qQa%'^ - ik'^q'^qoo^b + 2k'^q^qQa^ - ik"^ q"^ qoab"^ + ik'^q'^qoab - 2k'^q^qoa 
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^k'^q^qoh'^ - k'^q^qoh - Ak'^q^qoS^v'^z^{l - zf - 4:k^q^qoS^vz^{l - z){ux) 
-Ak'^q^qoSv'^z'^{l - z)[u{l - x)] - Ak'^q'^qoSvz'^{ux)[u{l - x)] + 2k'^q^a%'^m 
-2k'^q^a%m+k^q^a^rn-2k'^q'^ab^rn+2k^q'^abrn-k^q'^arn+k'^q%'^rn-k'^q%rn 
-2k'^q^ms'^v^z^{l-z)^-2k^q^ms^vz^{l-z){ux)-2k^q^msv^z'^{l-z)[u{l-x)] 
-2k'^q^msvz\ux)[u{l - x)] - APq^a^b^ + APq^a^b - 2k'^qla^ + APq^ab'^ 
-'ik'^qQub + 2k'^qla + 2k'^qla?m - 2k'^qlam - Ak'^qQO^w? + Ak'^qoaw? 
+2k'^qom^v^z^{l - z) + 2k^qom^v'^z{l - zf + 2k'^qom^v^z{l - z) 
+2k'^qQm^vz^{ux) + 2k^qQm'^vz(l — z)(ux) + 2k'^qQm^vz{ux) — k'^qom? 
-2k'^o?rrr' + 2k'^am^ + k'^m^v'^z'^{l - z) + k'^m^v^z{l - zf + A;^mV;2(l - z) 
-\-k^m^vz^{ux) + k^m^vz{l — z){ux) + k^m^vz{ux) — k^m^ 
-%{kqfqQs\'^z^{l - zf + l&{kqf q^sv z^ {I - z){ux)[u{l - x)\ 
-S{kqfqQZ^{uxf[u{l - x)f + A{kqfa%'^m - A{kqfa^hm - A{kqfab'^m 
+A{kqYabm — A{kqf'ms^v'^z'^{l — zY + S{kqYmsvz'^{l — z){ux)[u{l — x)] 
-A{kqfmz^{uxY[u{l - x)f + A{kq)q^qQab^ - 6{kq)q'^qoab'^ + 2{kq)q^qQab 
-2{kq)q''qob^ + 3{kq)q\ob'' - {kq)q^qob - 8{kq)q''qoS^vz\l - zf 
-S{kq)q^qoS^vz'^{l - z)[u{l - x)] + %{kq)q^qoS^ z^{l - z){ux)[u{l - x)] 
-\S{kq)(f q^sz"^ {ux)[u{\ — x)]^ + A{kq)q^ab'^m — Q{kq)q^ab'^m + 2{kq)q^ abm 
—2{kq)(f'b^m + 3{kq)q^b^m — {kq)q^bm — A{kq)q^ms^vz^{l — zf 
-A(kq)q'^ms^vz^{l - z)[u(l - x)] + A(kq)q^ms^z^(l - z){ux)[u(l - x)] 
+A(kq)q'^msz'^ {ux)[u{l — x)f — A{kq)qoabm'^ + 2{kq)qoam'^ + 2{kq)qobm'^ 
+A[kq)qo'm'^ sv z'^ {1 — z) + A{kq)qQm'^svz{l — zf + A{kq)qQm^svz{l — z) 
—A{kq)qo'm^ z^ {ux)[u{l — x)] — A{kq)qo'm^ z{l — z){ux)[u{l — x)] 
-A{kq)qom^ z{ux)[u{l - x)] - {kq)qom^ - A{kq)abm^ + 2{kq)am^ + 2{kq)bm^ 
+2{kq)m^svz^(l - z) + 2{kq)m^svz{l - zf + 2{kq)m^svz{l - z) 
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—2{kq)m^ z'^ {ux)[u{l — x)] — 2{kq)m^z{l — z){ux)[u{l — x)] 

-2{kq)m^z{ux)[u{l - x)] - {kq)m^ + 2g^go&^ - Vgo&^ + '2q%b'^ 

-2q%s^z\l - zf - 4g^goS^-s'(l - z) [u{l - x)] - 2g^gos'^'[ii(l - x)f + g^6^m 

-2q%^m + gVm - g^msS^(l - - 2q'^ms^z^{l - z)[u{l - x)] 

—q'^ms'^z^[u{l — x)]^ — Aq^qob'^m?' + Aq^qobm^ + 2(f' qQ-rn^ z^ [1 — z) 

+2q^qom^s^z{l — zY + 2q^qom^s^z{l — z) + 2q^qom^sz^[u{l — x)] 

+2q'^qom^sz{l-z)[u(l-x)]+2q^qorn^sz[u(l-x)]-2q'^qorn'^-2q%^rn^+2q%rn^ 

+g^m^s^2;^(l - z) + q^rrr's^z{l - zf + q^rrV's^z{l - z) + q^m^sz'^[u{l - x)] 

+q^rrV' sz{l — z)[u{l — x)] + (fm^szluil — x)] — (fm^ + 2qlm'^ — 2qlm^ . 

Next we go to calculation of the upper photon loop in Fig.lc. To facilitate 
transformations we introduce new notation which resembles the one used in 
the previous section 

= m^-z{k'^{l-ux)[l-z{l-ux)]+q^[l-u{l-x)][l-z{l-u{l-x))] (64) 
+2kq[l -u- z{l- ux){l - u{l - x))]] 
= + - 2/3^ (A;g) - ^,m\ 



where 



(1 — ux)[l — z{l — ux)] 



Iz = 



fi- 


u{l - 


-x)][l-z{l- 


u{l - 


^))]' 


ll 


— u — 


z{l — ux){l — 


u{l - 




[1- 


-u{l 


-x)\[l-z{l- 

1 


-u{l - 


-^))]' 


z[l- 


-u{l - 


-x)\[l-z{l- 


■u{l - 





Note that 



az{z — 1) — a, (65) 

Pz{z = 1)=P, 
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lz{z = 1) = 7. 

As was discussed in detail in the case of the ladder diagram we have to 
introduce small photon mass A for the photons in the upper loop to make all 
intermediate integrations infrared safe. 

Then 

(1 - t)[{l - y){q^ - A^) + y{q^ + 2mq,)] + t(-7.A,) = q"^ - rf^, (66) 
where 

q" = q + y{l-t)m- P,tk, 

dl, = -kH{a, - f3lt) + m2[y2(i _ tf + 7,t] + \\l -t){l-y). 
Hence, 

(g^-A^)^(,' + 2mgo)A. = "^^^ [ ^^^^ ' lo ' '^\,'- - dl^ 

(67) 

1 



2A,l[dyil-y)[dUil-tf' 



(g2 - A2)2(g2 + 2mgo)A2 7o ^'7o ' ' {q"^ - dl 



■\z) 



Next we perform euclidean rotation and calculate the integral over q'^ ^ 
for the contribution of the crossed diagram to the energy shift in eg. (p!0|) . 

AE, = 3 ^ — - hm 68 

vm'^ m TT"^ J Ank^ Jo 

//2 „//2 



/^rfwu ^ dz-f, [\y{l-y) dt{l-tf f dq"^q 
Jo Jo Jo Jo J 

2F^{q = q" -y{l-t)m + f5,tk) 2z^F^{q = q" - y{l - t)m + p,tk) 

Explicit expressions for the numerators (after Wick rotation) in the inte- 
grand have the form 

2F^ = -3q"'^mP^ + 2k'^mP^ + 2m^[y{l - t)fP^, (69) 

where 



^Subscript E is omitted below. 
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= 86'z2[?/(l - t)] - 2b'' - Sbz^[y{l - t)] + 2bz' (70) 
-65^/(1 - z)[y{l - t)] - 4sz'[u{l - x)Ml - t)] - 2s[u{l - x)Ml - t)] 
-2z'[u{l - x)ny{l - t)] + 4z'[u{l - x)][y{l - t)] - 2z'[y{l - t)] 
+2/(1 - z)[u{l - xmi - t)] - Az^l - z)[u{l - x)][y{l - t)] 
+2/(1 - zMl - t)] + z'[y{l - t)] - / + 2z[uil - x)ny{l - t)] 
-4z[u{l - x)Ml - t)] + 2z[y{l - t)] - 2(1 - z)[u{l - x)]'[y{l - t)] 
+4(1 - z)[u{l - x)][y{l - t)] - 2(1 - zMl - t)], 

= 8a2/[|/(l - t)] - 2a^z^ + I2ahz'[y{l - t)](/?,t) (71) 
-iabz'^iP^t) - 6a/[y(l - t)]{p,t) - 8a/[y(l - t)] + 2az'{(5^t) + 2a/ 
+10&2/[y(l-i)](/3,i)2-262/(/3,i)2-106/[y(l-i)](/3,i)2-6&/[y(l-i)](/3,i) 
+2bz\^tf + 2bz\f3,t) - 6//(l - z)[y{l - t)]{f3,tf 
-I2svz\l - z)[y{l - t)]{P,t) - Asz^[u{l - x)][y{l - t)]{l3,tf 
-2s[u{l - x)Ml - t)]{(3,tf - Qv\\l - z)[y{l - t)] - Avz\ux)[y{l - t)] 
-2v{ux)[y{l - t)] - 2z\uxf[y{l - t)] - Az\ux)[u{l - x)][y(l - t)]{l3,t) 
-rAz\ux)[y{l - t)]{^t) + Az\ux)[y{l - t)] - 2z'[u{l - x)f[y{l ^ 
+4/[«(l-a;)]b(l-t)](/5,t)2 + 4/[t.(l-x)][y(l-t)](/3,t)-2/[y(l-t)](/3,t)2 
-4/[t/(l - t)]{(5,t) - 2/[y(l - t)] + 2/(1 - z){uxf[y{l - t)] 
+4/(1 - z){ux)[uil - x)Ml - t)]{f3,t) - 4/(1 - z)iuxMl - t)]{f3,t) 

-4/(1 - z){uxMl - t)] + 2/(1 - z)[u{l - x)ny{l - t)]{P,tf 
-4/(1 - z)[u{l - x)Ml - mp^tf - 4/(1 - zMl - x)][y{l - t)]{P,t) 
+2/(1 - z)[y{l - t)]{P,tf + 4/(1 - z)[y{l - t)]{(3,t) + 2/(1 - z)[y{l - t)] 
+8z\ux)[u{l - x)Ml - t)]{P,t) + 2/[y(l - t)]{P,tf + 3/[y(l - t)]{P,t) 
+/[y(l - t)] - z\(5,tf - z\(5,t) - / + 2z{uxf\y{\ - t)\ 
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+Az{ux)[u{l - x)Ml - - Az{ux)[y{l - t)](/3,t) - Az{ux)[y{l - t)] 

+2z\u{l - x)f[y{l - t)mtf - Az[u{l - x)Ml - t)mtf 
-4z[M(l-x)][y(l-t)](/3,t) + 2z[y(l-t)](/5,t)2 + 4^[y(l-t)](/?,t)+2^[l/(l-t)] 
-2(1 - z){uxf[y{l - t)] - 4(1 - z){ux)[u{l - x)][y(l - t)](/5,t) 
+4(1 - z){ux)[y{l - t)](/3,t) + 4(1 - z){ux)[y{l - t)] 
-2(1 - z)[u{l - x)f[y{l - t)]{(3.tf + 4(1 - z)[u{l - x)][y(l - t)mtf 
+4(l-z)[n(l-a;)][y(l-t)](/5,t)-2(l-^)[y(l-t)](/5,t)2-4(l-z)[y(l-t)](/3,t) 
-2(1 - z)[y{l - t)] + A{ux)[u{l - x)Ml - t)]{(3,t), 

-Gs'^z^il - z)[y{l - t)] - Asz^[u{l - x)][y{l - t)] - 2s[u{l - x)][y{l - t)] 

-2z'[u{l - x)ny{l - t)] + 4z'[u{l - x)Ml - t)] - 2z'[y{l - t)] 

+2^3(i_^)[^(l_^)]2[^(^_^)]_4^3(^_^)[^(^_^)][^(^_^)]_^2z3(l-^)[y(l-t)] 

-3^2 + 2z[u{l - x)f[y{l - t)] - Az[u{l - x)] [y{l - t)] + 2z[y{l - t)] 

-2{l~z)[u{l-x)f[y{l-t)]+A{l-z)[u{l-x)Ml-t)]-2{l-z)[y{l~t)l 

Second term in the integrand in eq.(p^) has the form 

2FI = 2q"^mTl - q"^PmT^ - q"^m^T^ (73) 

+2k^mT^ + 2k'^m^T^ + 2m^[y{l - t)]'^T^, 

where [] 

= -4b'^[y{l - t)] +b^ + 8b^[y{l - t)] - 2b^ (74) 
-46'[l/(l - t)] +b^ + As*z\l - zflyil - t)] - sV(l - z)^ 
+8sh'^{l - z)[u{l - x)][y{l - t)] - 2sV(l - z)[u{l - x)] 
+4sV[u(l - x)]^[y{l - t)] - s^z^[u{l - 



^It is not difficult to check that despite its cumbersome appearance function vanishes 
identicaUy: r| = 0. 
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= -6a%'^[y{l - t)] + 6a%^ + 6a%[yil - t)] - 6a% (75) 

-3a^[y{l-t)] + 3a^ - 16ab%{l - t)]{p,t) + 12ab^{p,t) +2Aab^[y{l - t)]{p,t) 

+6ab'^[y{l - t)] - 18ab'^{p^t) - 6ab'^ - 8ab[y{l - t)]{p,t) - 6ab[y{l - t)] 

+6ab((3,t) + 6ab + 3a[y(l - t)] - 3a - 166^[?/(1 - t)]((3,t)^ + 6b\(3,t)'^ 

+32b'[yil - t)]{(3,tY + 863[y(l - t)]{P,t) - Ub^p^tf - 6b%P,t) 
-1662[y(l-t)](/3,t)2-i262[^(l_t)](/3,t)-362[^(l-t)] + 662(/5^t)2 + 952^^^^) 

+262+46[y(l-i)](/3,i)+36[y(l-i)]-36(/?,t)-26+16sS2(l-^)2[y(l-i)](/?,i)2 

-6sS'(l - zfiP.tf + 32sW(l - z)^[y{l - t)]{P,t) - 12sW(l - zfiP.t) 

+32s'z'il - z)[u{l - x)Ml - t)]{P,tY - Us'z'il - z)[uil - x)]{P,tY 
+16s\^z^{l - zf[y{l - t)] - 6s\2^^(l - zf + 12s^vz^{l - z){ux)[y{l - t)] 

-4sW(l - z){ux) + 32sW(l - z)[u{l - x)][y{l - t)]{p^t) 
-12s\z'^(l - z)[u(l - x)](P,t) - 32s^z^(l - z)(ux)[u(l - x)][y(l - t)]{P,t) 
+12sV(l - z){ux)[u{l - x)]{(3,t) + 16sV['u(l - x)f[y{l - t)]{(3^tf 
-6s'^z^[u{l - x)]'^{P,tf + 12st;V(l - z)[u{l - x)][y{l - t)] 
-^sv^z^{\ - z){u{\ - x)\ - 8svz^{\ - z){ux){u{\ - x)\{y[\ - t)\ 
^^svz^{\ - z){ux)\u{\ - x)\ + \2svz^{ux)\ii{\ - x)][y(l - t)\ 
-\svz^[ux)\u{\ - x)\ - 32sz^{ux)\u(l - x)Y\y{\ - t)]{p,t) 
+12sz'^{ux)[u{l - x)]'^{p^t) + Az^{uxY[u{l - x)f[y{l - t)] 
-2z\uxf[u{l - x)]^, 

= -m^[y{l - t)f + Qb'[y{l - t)f + 3263[^(i - t)f (76) 
-I2b^[y{l - t)f - m''[y{l - t)f + 662[y(l - t)f + I2b''[y{l - t)] - Ab^ 
-12fe[y(l - t)] +4b+ 16sS'(l - z)^[y{l - t)]^ - 6sV(l - ^)'[?/(l - t)]^ 
+32sV(l - z)[u{l - x)][y{l - t)f - 12sV(l - z)[u{l - x)][y{l - t)]^ 
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-6s^z^{l - z)[y{l - t)] + 2s^z'^{l - z) + 16s^z^[u{l - x)]^[y{l - t)]^ 
-6s^z^[u{l - x)]^[y{l - t)]^ - 6s^z{l - zf[y{l - t)] + 2s^z{l - zf 
-6s^^(l - z) [y{l - 1)] + 2s^z(l -z)- Qsz^ [u{l - x)\[y{l - 1)] + 2sz\u{l - x)] 
-6s^(l - z) [u{l - x)] [y{l - 1)] + 2sz{l - z) [u{l - x)] - Qsz[u{l - x)] [y{l - 1)] 

+2sz[u{l - x)] + 3[2/(l - t)] - 3, 

= -2a%{l - t)] + - Aa%y{l - t)]if3,t) + Aa%{l3,t) (77) 
+2o^[y{l - t)]{l3,t) + Aa'[y{l t)] 2a'{f3,t) - 2a' - Aa%''[y{l - t)]{l3,tf 
+6aV(/5^t)2+4a26[y(l-t)](/5^t)^ + 6a^6[l/(l-t)](/3,t)-6a^6(/3,t)=^-6a26(/3^^ 
-2a^[y{l - t)]{P,tf - 3a'[y{l - t)]{/3,t) - 2a'[y{l - t)] + a\P,tf + ?,a\P,t) 
W-Aah^[y{l-tW,tf+Aah\p,tf+^ah-'[y{l-t^^^^^ 

-6ab\P,tf - Gab^P^tf - 2ab[y{l - t)]{P,tf - 4a6[y(l - t)]iP,tf 
-2ab[y{l - t)](,4t) + 2ab{l3,tf + Qab{l3,tf + 2ab{(3,t) + 2a[y{l - t)]{(3,tf 
+a[y(l - tW^t) - aiP^tf - a{P,t) - 2b^[y{l - t)]{(3,t)' + b^P^t)' 
+Ab'[y{l - t)mtr + 2b'[y{l - t)mtf - 2b\(5^tY - 2b\(5^tf 
-2b''[y{l-t)]{^tf-?,b''[y{l-t)]{P,t)'-b^^^ 

+b\p,tf + b[y{l - mp^tf + b[y{l - mp^tf - b{P,tf - b{P,tf 
+2sV(l-^)2[y(l-i)](/?,i)^-sV(l-^)2(/?,i)^+8sW(l-^)2[y(l-i)](/?,i)^ 
-4sW(l - zfiP^tf + As^z^il - z)[u{l - x)Ml - t)]{P,tY 
-2s^z\l - z)[u{l - x)]{P.tY + 12sVz^(l - zf[y{l - t)]{P,tf 
-Q.s'v^z^il - zf{p,tf + As\z^{l - z){ux)[y{l - t)]{P,tf 
-2sW(l - z){ux){p,tf + 8sW(l - z)[u{l - x)][y{l - t)]{p,tf 
-4sW(l - z)[u(l - x)](P^tf - 8s^^^(l - z)(ux)[u(l - x)][y{l - t)](P,tf 

+As^z\l - z){ux)[u{l - x)]iP,tf + 2s''z''[u{l - x)]^[yil t)]{PJ)^ 
-s'z'[u{l - x)]\p,ty + 8sv'z\l - zf[y{l - t)]{P.t) - isvh'{l - zf{P,t) 
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+8st;V(l - z){ux)[y{l - t)]{P,t) - U'v^z^iX - z){ux){p,t) 
+isvh''{l - z)[u{l - x)][y{l - t)]{P,ty - 2sv'z^{\ - z){u{\ - x)\{^,tf 
-\^svz^{\ - z) (ux) [u{l- x)] [y{l- t)]{P^tY + 8svz'^{l - z) {ux) [u{l - x)\{(5^tf 
+Asvz'^{ux)[u{l - x)][y{l - t)](/3^t)^ - 2svz^{ux)[u{l - x)\{(3^tf 
-Ssz^{ux)[u{l - x)]'^[yil - t)]{p,tf + Asz^{ux)[u{l - x)f{p,ty 
+2v^z\l - zfiyil - t)] - v^z^l - zf + ^v''z^{\ - z){ux){y{l - t)\ 

-2v^z^{l - z){ux) - Sv'^z'^il - z){ux)[u{l - x)][y{l - t)]{P^t) 
+iv^z'^{l - z){ux)[u{l - x)]{p,t) + 2v'^z'^{uxf[y{l - t)] - v^z^{uxf 
-Svz^{uxf[u{l - x)][y{l - t)\{(3,t) + Avz^{uxf[u{l - x)]{(3,t) 
-rSz^uxfH-^ - ^)fW - t)mtf - Az^uxfH-^ - x)f{(5S\ 

n = 2a%''[y{l - t)f - 2a%[y{l ~ f)f + ^a^[y{l - t)f (78) 
+4a2[|/(l - t)] - 2^ - Aa}?{y{\ - t)\\i5,t) + Aa}?{y{\ - t)\\(i,t) 
+6a62[y(l-i)]=^(/3,i)-6a62[y(l-i)]^(/3,i)-2a62[y(l-i)]2-2a6[y(l-i)]=^(/3,i) 
+2a%(l-i)]2(/3,i)+2a6[y(l-i)]2+4a%(l-i)](/3,i)-4a6(/3,i)-3a[z/(l-i)]2 
-2a[i/(l - t)](/5,t) - 4a[y(l - t)\ + 2a(/3,t) + 2a - %%{\ - t)f{^,tf 
^2h%[\ + 86^^(1 - t)\\{i-.tf + 2\?\y{\ - tyfm) 

-Ah\y{\ - t)np,tf - 2b'[y{l - t)rm - Ab'[y{l - t)]'m' 
-3b'[y{l - t)f{P,t) - b'[y{l - t)r + 2b'[y{l - t)r{P,t)' + 3b'[y{l - 
+62[^(l-t)]2 + 462[^(l_t)](;5,t)^-262(/3,t)2 + %(l-t)]3(/3,t)+6[y(l-t)]3 
-6b(l-t)]2(/5,t)-6[l/(l-t)]2-46[y(l-t)](/?,t)2-26[y(l-t)](/5,t) + 26(/5,t)2 
+2b{P,t) + 4sV(i _ ^y^y^i _ t)f{p,tf - 2s'z\l - zf[y{l - t)f{(5,tf 
+%s\z\l - zf[y{l - t)fm - As\z\l - zf[y{l - t)fm 
+8sV(l-^)[w(l-x)][|/(l-t)]3(/5,t)2-4sV(l-;.)[«(l-x)][|/(l-t)]2(/3,t)2 
+As\^z\l - zf[y{l - t)f - 2sV^2(l - zf[y{l - t)f 
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+4sW(l - z){ux)[y{l - t)f - 2sW(l - z){ux)[y{l - t)]^ 
+8sW(l-^)[K(l-x)][y(l-i)]3(/?,i)-4sW(l-^)[K(l-x)][y(l-i)]2(/?,i) 
-Ss^z'il - z)iux)[uil - x)Ml - t)]^iP,t) 

+4sV(i - z){uxMi - x)Mi - t)]\p,t) - 2sV(i - zMi - t)]{p,ty 

+sV(l - z){(3,tf + 4s'z'[u{l - x)ny{l - t)]'{m' 
-2sV[«(l-x)]2[y(l-t)]2(/?,t)2-2s2^(l-;.)2[y(l-t)](/3,t)2+s2^(l-^)2(/3,t)2 

-2s2z(l-z)[y(l-t)](/3,i)2+s2^(l-z)(/3,t)2+4s^;V(l-z)[ii(l-x)][y(l-i)]^ 
-2sv^z^{l-z)[u{l-x)][y{l-t)]'^-Asvz^(l-z)[y(l-t)]{P,t)+2svz^{l-z)(P:,t) 

+Asvz'^{ux)[u{l - x)][y{l - t)f - 2svz'^{ux)[u{l - x)][y{l - t)f 
-4svz{l - zfiyil t)]{0,t) + 2svz{l - zfiP.t) - Asvz{l - zMl - t)]i(3,t) 
+25^^(1 - z){p;,t) - 8sz'^{ux)[u{l - x)f[y{l - t)]^{Pzt) 
+Asz\ux)[u{l - x)ny{l - t)r{P,t) - 2sz'[u{l - x)Ml - mp.tf 
+sz^[uil - x)]i(3,ty - 2sz{l - z)[uil - x)Ml - t)]if3,ty 
+sz{l - z)[u{l - x)]iP,tf - 2sz[uil - x)][yil - t)]{P,tf 
+sz[u{l - x)]{p,tf - 2v^z\l - z)[y{l - t)] + v^z^l - z) 
-2v^z{l - zflyil - t)] + v'^z{l - zf - 2v^z{l - z)[y{l - t)] + v^z{l - z) 
-2vz^{ux)[y{l - t)] + vz'^{ux) - 2vz(l - z)(ux)[y{l - t)] + vz{l - z){ux) 
-2vz{ux)[y{l -t)] +vz{ux) + Az"^ {ux)[u{l - x)][y{l -t)]{P^t) 
-2z'^{ux)[u{l - x)]{p,t) + 4^(1 - z){ux)[u{l - x)][y{l - t)]{p,t) 
-2z{l - z){uxMl - x)]{f3,t) + 4z{ux)[u{l - x)][y{l - t)]{P,t) 
-2ziux)[u{l - x)mt) + 2[y{l - + [y(l - + [y(l - t)] 

- iPzt) - 1, 

= -2b'[y{l - t)f + b'[y{l - t)]' + 4b'[y{l - t)f (79) 
-2b'[y{l-t)]'-2b'[y{l-t)]' + b'[y{l-t)]' + 4b'[y{l-t)]-2b'-4b[y{l-t)] 
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+26 + 2sV(l - zf[y{l - t)f - s^z\l - zf[y{l - t)]^ 
+4sV(l - z)[u{l - x)][y{l - t)f - 2sh''{l - z)[u{l - x)][y{l - t)]'' 
-2s2^2(l - z)[y{l - t)] + s^z^il -z)+ 2s^z^[uil - - t)f 

-sV[m(1 - x)]^[y{l - - 2s^z{l - zf[y{l ~ t)] + s^z{l - zf 
-2.s^z{l - z)[y{l - t)] + s'^z{l -z)- 2sz'^[u{l - x)\[y{l - t)] + sz'^[u{l - x)] 
-2sz{l - z)[u{l - x)] [y{l - t)] + sz{l - z)[u{l - x)] - 2sz[u{l - x)] [y{l - t)] 

-\-sz[u{l — x)] — 3. 

After integration we obtain 

a''{Zaf 



AEc = — ^-^mi—f— / dx / duu / dz-^^ (80) 
Trn"* m Jo Jo Jo 

J\y{l-y)l'dt{l-trj ^{47.tk^ 



^mz^T^ 



3mP[= m^[y{l-t)fP^ mz'^Tl 



3rfL " 2dl 

^ m^z^T^ , m^z^\y(l-t)fT;^, 

-2-f,t -r^ + A-i,t — '-^—^}. 

Mi, Ml, ^ 

Integration over |k| is no more difficult than in the ladder case. We 
introduce notation which resembles the one used in the ladder case 

dl, = ViHia, - Pit) + m^[y\l - tf + ^,t] (81) 

^ p^^^2 ^ m'[y\l-tf + ^^t] + \\l-t){l-y) , 

Pz 



where 



p, = t{a,-(3lt), (82) 
m^[y2^l-tf^^A-r\\l-t){l-y) 



\z 
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Once again as in the ladder case we encounter the problem of apparent 
divergence of the integral for the energy shift. As was already explained it 
is necessary to check that the total electron-line factor vanishes as k^ in the 
low frequency region. It is an easy task to do. We simply omit integration 
over k in eq. (^) and eq. (|80|) and put k to be equal to zero in the integrand. 
Then we obtain 

Aetest = Aei + Ae^ (83) 



1 rl rl rl 

dx du dy{l - y) / dt{l - tf 
Jo Jo Jo 



{l-u)Pl -t)]2(l -m)PJ -m)T, 



[u{l - x)]u{l - u)Ti [u{l-x)]u[y{l-t)f{l-u)Ti ^ 



^Jo "^'^'^ 3[y2(l-t)2+^^t] +3[y2(l_t)2+^^i]^ 



2[y2(i_t)2 + ^^t] - 6[y2(l-t)2 + ^^t]2 
z''u[y{l-t)]^T^ 

We calculated this integral numerically and it turned out to be equal 
zero as expected. Then there are really no difficulties in the treatment of 
the integrals contained in eq. (|80|) ; one has simply to perform subtraction of 
the leading low frequency terms in the electron-line factor since these terms 
cancel in any case in the total electron-line factor. 

Subtraction of the low frequency part is necessary in the last five terms 
of the integrand in eq.(^). The numerators of these terms are independent 
of momentum k and subtraction may be easily performed with the help of 
the identities 

^ 1_ _ k^ ^ k^ 

11 k2 k2 
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11 k2 k2 



k2 k2 k2 



2 



Substituting these subtractions in eq.(|80D and performing the momentum 
integration we obtain 

Ai^r^ = "'^^f f'dx (85) 

irn-^ m n Jo 



1 rl rl rl 

duu I dz'jz / dy{l — y) dt(l — t)' 
Jq Jo Jo 



Izt ,^^,2 .^c .2^c m^z^[y{l-t)fT^ 



r '2'' r 2 2t-ic 3 2t-ic 

^Pz^Az ^Az 

1 ^oc,,,.^.2^c m3[y(l-t)]2p| 



^^Pi^xz Pz 



\z 

6 1 



+izt — - 27,t — 4 ^ 

^Az Pz^lz 

+7 ^ 3Pi^ 2 - 67.^^ 

6p.^A^ P.^A^ 

5 TOTAL CONTRIBUTION TO THE LAMB 
SHIFT AND DISCUSSION OF RESULTS 

It is straightforward to calculate numerically the total contribution to the 
S'-level Lamb shift of order a'^{Za)^m induced by the light by light scat- 
tering insertion. This contribution is given by the sum of the expressions 
in eq.(|5UD and in eg. (|53D . The only subtlety is that these expressions con- 
tain an auxiliary infrared regularizing parameter A. As we have mentioned 
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above the final formula for the energy shift should be given in the limit of 
the vanishing photon mass and should produce an unambiguous value for the 
contribution to the Lamb shift. We have checked numerically that the result 
of integration is independent of this small intermediate photon mass as it is 
varied from = 10~^ to = 10~^. This limited range of values for the 
auxiliary photon mass is defined by the hmited accuracy and productivity 
of the computer and may be widened at the expense of significantly more 
computer time. Our final result is 

a^{Za)^ ,m. 



AE = -0.122(2) 



r \3 



m. 



m 



(86) 



One may readily obtain more digits in the expression above if necessary. 

In this work, initiated in |]l| , we have finished calculation of all corrections 
to the Lamb shift of order a'^{Za)^m induced by the diagrams with closed 
electron loops. Our results are presented in the Table.0. 



Corrections to the Lamb shift 
induced by the diagrams with closed electron loops 



AE 




2^ (kHz) 


IS (kHz) 


afl 


-0.061 


-0.329 


-2.63 


fei 


0.508 


2.747 


21.98 


c§ 


0.611 


3.305 


26.44 


d§ 


-0.073 


-0.394 


-3.15 


e-this work 


-0.122(2) 


-0.660 


-5.28 


Total 


0.863(2) 


4.67(1) 


37.3(1) 



The theoretical contributions to the Lamb shift presented in the Table 
are clearly necessary for comparison of theory with recent experiments on 
measurement of the 25*1/2 — splitting in hydrogen ||^, |^ 

Au = 1057 845 (9) kHz , 

®While this paper was in preparation a wor k |p[ appeared which contains recalculation 
and confirmation of the results obtained in ||l|, ^ ^ as well as the first calculation of the 
light by light contribution in a framework completely different from the one used in this 
paper and in || . The result for the light by light contribution presented above is in 
agreement with the respective result in 
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Az/ = 1057 851.4 (1.9) kHz, 



(87) 



and on measurement of the IS* Lamb shift ^ 

Az/15 = 8 172.82 (11) MHz. (88) 

Work on the calculation of the still unknown corrections of order a'^{Za)^m, 
produced by the diagrams with two radiative photon insertions in the elec- 
tron line is in the finishing stage now by at least two groups|. 
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Figure Caption 



Fig.l. Gauge invariant set of diagrams with light by hght scattering insertions 
in the external Coulomb lines. 
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